Condensate fraction and atomic kinetic energy of liquid 3 He- 4 He mixtures 
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We present neutron scattering measurements of the momentum distribution of liquid 3 He- 4 He 
mixtures. The experiments were performed at wavevectors Q, 26 < Q < 29 A -1 , on the MARI time- 
of-fiight spectrometer at the ISIS Facility, Rutherford Appleton Laboratory, a spallation neutron 
source. Mixtures with 3 He concentrations x between and 20% were investigated both in the 
superfluid and normal phases. From the data, we extract the Bose-Einstein condensate fraction no 
and the momentum distributions of 3 He and 4 He atoms. We find that no increases somewhat above 
the pure 4 He value when 3 He is added; e.g from no = (7.25 ± 0.75)% at x — to (11 ± 3)% at x = 
15-20%. This agrees with predictions but is less than the only previous measurement. We find a 4 He 
kinetic energy K4 for pure 4 He that agrees with previous determinations. K4 decreases somewhat 
with increasing 3 He concentration, less than observed previously and found in early calculations but 
in agreement with a more recent Monte Carlo calculation. The 3 He response is not well reproduced 
by a Fermi gas momentum distribution, n(k). Rather an n(k) having a small step height at the 
Fermi surface and a substantial high momentum tail characteristic of a strongly interacting Fermi 
liquid provides a good fit. This n(k) is consistent with calculated n(k). Thus agreement between 
theory and experiment is obtained comparing n(k) in contrast to earlier findings based on comparing 
calculated and observed 3 He kinetic energies. 

PACS numbers: 67.60.-g, 61.12.Ex 



I. INTRODUCTION 

Liquid 3 Hc- 4 He mixtures are excellent examples of in- 
teracting Bose and Fermi liquids in nature. Their behav- 
ior is dominated by quantum-mechanical exchange effects 
which involve not only identical particles, as in the pure 
liquids, but also interactions between particles obeying 
different statistics. As a result, a rich array of macro- 
scopic properties are observed as the 3 He concentration 
is increased. Among these properties are a gradual sup- 
pression of the 4 He superfluidity and a finite solubility 
(~6.6%) of 3 He in 4 He at absolute zero temperature. 

Historically, investigationsi*2i2i4 of this binary mixture 
involved thermodynamical and hydrodynamical experi- 
ments. A central goal was to measure the phase diagram 
of the mixture and to investigate the elementary excita- 
tions as a function of 3 He concentration, x. In particular, 
the 3 He quasiparticles and the 4 He phonon-roton (p-r) 
spectral behavior have drawn much attention££iL& 

Fifty years ago, Landau and Pomeranchuk (L-P) pro- 
posed a model of the dilute helium mixtures at low tem- 
peratures as two interpenetrating fluids of 3 He and He 
atoms for which 3 He- 4 He interactions are retained and 
3 He- 3 He interactions are ignored. In this simple picture, 
the 3 He quasiparticle energy is given by es(k) = Hk 2 /2m* 
where m* is the 3 He effective mass in the presence of 4 He 
atoms. At low k, the 3 He excitation spectrum observed 
in neutron scattering^ agrees well with the L-P spectrum. 
However, for k > 1.5 A -1 , it departs from the theory and 
falls below the L-P curve. This observation and the ob- 



served small shift in the 4 He roton energji^, A, suggest a 
dynamical level-repulsion between es(k) and 64 (fc) in the 
vicinity of A. A smaller shift in the roton energy was also 
observed in Raman scattering 7 - for several mixtures con- 
centrations, up to x = 30%. As a possible explanation, 
Pitaevski: 9 proposed the existence of a small roton-like 
minimum in €3 (k). This suggestion found some theoreti- 
cal support^ but no experimental evidence. 

A fundamental parameter of interest in helium mix- 
tures is the energy-dependent quasiparticle effective 
mass, m*. Yorozu and collaborators^ reported high ac- 
curacy measurements of the 3 He effective mass in mix- 
tures and found an m* that strongly depends on pressure 
but varies little with x. Their result was later supported 
by calculations^ Recently, the effect of m* on the dy- 
namics of pure liquid 3 He has been investigated^ 

The first neutron scattering measurements at low mo- 
mentum transfer, Q, aimed at probing elementary exci- 
tations in helium mixtures were reported by Rowe et al& 
for a 5% mixture at saturated vapor pressure (SVP) and 
a temperature of 1.6 K and by Hilton et al& for several 
mixtures at SVP for temperatures in the range 0.6-1.5 K. 
In these pioniering investigations, the cryogenic capabil- 
ities were limited to temperatures close to, or above, the 
Fermi temperature, Tp . In this limit, thermal broaden- 
ing is important and can mask the quasiparticles energy 
band. Both Rowe et al. and Hilton et al. reported a 
small shift of the 4 He p-r curve but disagreed on the sign 
of the shift. 

Fak et al& were the first to investigate the excita- 
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tions at temperatures below Tp. Their measurement 
confirmed the small shift in the p-r spectrum observed 
by Rowe et al. and by Hilton et at. They found no 
evidence of a roton-like minimum in the 3 He spectrum. 
Their experimental results are consistent with the theo- 
retical predictions of Bhatl*^ and of Gotze et al.w^ 

Other complementary experimental techniques on he- 
lium mixtures include several studies of the elementary 
excitations by Raman scattering 7 ^ and a study of the 
static structure factor S(Q) by X-ray scattering^. We re- 
fer to the review articles by Ahlers 3 , Baym and Pethick 18 
and Glyde and Svensson-S and to the recent book by 
Dobbs 20 for further details on the subject. 

In recent years, inelastic neutron scattering at high 
Q has been widely used to measure the conden- 
sate fraction and atomic kinetic energy of pure liquid 
4 Hei 21 ' 22 i 23 i 24 i 25 ' 26 These high Q measurements are very 
challenging because of limitations imposed by Instrumen- 
tal Resolution (IR) and Final-State interactions (FS). A 
general review of the method and experiments on the 
VESUVIO neutron spectrometer is presented by An- 
dreani et al.SL In experiments involving 3 He, the large 
3 He neutron absorption cross-section reduces further the 
scattering intensity and thus the statistical precision 
of the experimental data. Nevertheless, a number of 
measurements at high Q have now been performed on 
pure liquid 2S i 29 i 30 i 31 ' 32 and solid 3 He^ 3 - In these measure- 
ments, sample cells having advantageous geometries with 
thin samples or samples at very low 3 He concentrations 
are generally used to minimize the effect of absorption. 

In contrast to the pure liquids, there are few reported 
neutron scattering measurements at high Q of the iso- 
topic helium mixtures. Most report values for the single- 
particle kinetic energies as a function of 3 He concentra- 
tion, x, in or close to, the normal phase 2 ^ 3 ^^ 1 .. The only 
neutron scattering measurement of no, to date, was re- 
ported by Wang and SokolM in a 10% mixture where 
they found an no = 18% ± 3. Several calculations of 
the Bose-Einstein condensate fraction (BEC), jiq, in he- 
lium mixtures 3 ^* 3 ^ 3 ^ 3 ^ 9 . predicted an enhancement of 
Uq. However, the enhancement predicted by these calcu- 
lations is significantly less than reported by Wang and 
Sokol. 

Theoretical treatments predict that the kinetic 
energy of both isotopes decreases with increasing 
concentrationi 35 i 39 i 40 i 41 Neutron scattering result a 29 ' 30 ! 31 
agree both quantitatively and qualitatively with theory 
for the behavior of if 4. In pure liquid 3 He, there is a 
reasonable agreement between recent experiments 29 ' 30 ' 31 - 
and theory 4 ^ on the value of K3, although earlier ob- 
served values were lower than calculated values^ 3 . Maz- 
zanti et a/J-L have recently obtained good agreement be- 
tween theory and experiment by comparing a calculated 
dynamic structure factor J(Q,y) directly with the ob- 
served J(Q, y). 

In dilute helium mixtures, in contrast, there is a sig- 
nificant disagreement in K3. For example, the calculated 
value of K3 for a 6.6% mixture is about 19 K [Ref. [3|| but 



the measured value (^10-12 K) is found to be the same 
as in the pure liquid 3 He independent of x. Repeated re- 
finements of the calculation technique using methods as 
diverse as Monte Carlo diffusion calculations^ 5 , to path 
integral MC techniques^ have failed to resolve this dis- 
agreement. 

In this paper, we report on deep inelastic neutron 
scattering measurements aimed at determining no, if 4, 
and the 3 He momentum distribution in 3 He- 4 He mix- 
tures. The quantity observed is the dynamic structure 
factor (DSF) S(Q, ui). For a sufficiently large momentum 
transfer from the neutron to the struck atom, Q — > 00, 
the observed DSF reduces to the impulse approximation 
(IA) 2 ^, 

S IA (Q, u) = j dkn(k)S(cu -us.- k.w R ). (1) 

where oj r — HQ 2 /2m and v R = HQ/m are the free atom 
recoil frequency and velocity, respectively. In the IA, 
Sia(Q,w) depends only on a single 'y scaling' variable 
y = (u> — ll>r)/vr and is conveniently expressed as 



Jia(v) = v R Si A (Q,u) = J dkn(k)S(y - k Q ) 



(2) 



where fcg = k.j^r. JjA{y) is denoted the longitudinal 
momentum distribution. At finite Q, however, the struck 
atom does not recoil freely but rather interacts with its 
neighbors. These interactions introduce a Final-State 
(FS) b roadening function R(Q,y) in J(Q,y) (see, e.g. 

for a detailed account of FS effects). Including 
these final-state interactions (FS), the exact J(Q,y) at 
high Q is 



J(Q,y)= / dse- isv J(Q, s) = / dse iys JiA(s)R(Q,s) 

. (3) 

where J(Q, s) is the intermediate scattering function, 
Jia(s) is the Fourier transform (FT) of Jia{v) and is 
the one-body density matrix (OBDM) for displacements 
along Q, and R(Q, s) is the FT of the FS function 
R(Q,y). The reader is referred to Boronat et al~ and 
Mazzanti et alM> for a detailed discussion of S(Q,u>) for 
the isotopic helium mixtures. 

When there is a condensate, the fraction of atoms with 
zero momentum transfer, no, contributes an unbroadened 
peak to Jia(jj) and can be directly extracted from the 
data. A central goal is to find an appropriate model for 
the momentum distribution n(k) including a condensate. 
For quantum systems, n(k) is also not a Gaussian. In 
earlier work 2 ^ 5 ^, we found that n(k) in liquid He and 
in solid 4 He differs markedly from a Gaussian. 

This report is laid out as follows. The experimental 
method is presented in section[Hj followed by a discussion 
on the data reduction and analysis methods in section 
IIIII The results are presented in section llVl followed by 
a discussion and conclusion in section Ivl 
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FIG. 1: MARI Data showing J(Q, y) for pure liquid 4 He in the nor- 
mal (open circles) and superfluid phase (closed circles). A higher 
peak arising from the Bose-Einstein condensate fraction is observed 
in the superfluid phase. The instrument resolution function is over- 
plotted (dotted line) for comparison. 



A special sample cell having a slab geometry, of the 
type described by Sokol et al£& and by Fak et al&, was 
designed. The sample cell was then placed in the beam 
and cooled using a 3 He sorption cryostat. Two separate 
experiments were performed on MARI under identical ex- 
perimental conditions. For each set of measurements, a 
separate background measurement of the empty cell was 
taken, and the two were found to be very similar, as ex- 
pected. Data was collected at temperatures of T = 0.4 K, 
T = 1.3 K and T = 2.5 K which were measured using Ge 
temperature sensors located at the bottom and top of the 
sample cell and connected to a Neocera temperature con- 
troller. A gas handling system at room temperature was 
used to prepare mixtures of 3 He gas to yield liquid 3 He 
concentrations of x — 0, 5, 10, 15, and 20 %. The pres- 
sure above the liquid was maintained at saturated vapour 
pressure, and was monitored with a pressure transducer 
by keeping a capillary line open between the gas handling 
system and the cell. 



B. Data Reduction 



II. EXPERIMENTAL DETAILS 

A. Experiment 

The experiment was performed using the MARI time- 
of-flight (TOF) spectrometer at the ISIS spallation neu- 
tron source at the Rutherford Appleton Laboratory, 
United Kingdom. MARI is a direct geometry chopper 
spectrometer in which the time of arrival of a neutron 
in the detector, measured from when the neutrons leave 
the moderator, determines its energy loss or gain after 
scattering from the sample. The momentum transfer de- 
pends on both the TOF of the neutron and its scattering 
angle. More than 900 3 He gas detectors provide a cover- 
age of scattering angles between 3° and 135° in steps of 
0.43°. A large range of momentum and energy transfer 
can therefore be observed simultaneously. 

An incident neutron energy of 765 meV was selected 
allowing wavevector transfers up to Q = 30 A" 1 and en- 
ergy transfers up to E = 700 meV to be investigated. 
Measurements of J(Q,y) were made for wavevectors 20 
< Q < 29 A -1 . The low energy resolution Fermi chop- 
per was selected to increase count rate in the presence 
of the highly absorbing 3 He. This has an energy resolu- 
tion of approximately 25 meV. The resolution function 
(in y) was determined at each Q as follows. The intrin- 
sic J(Q,y) of pure liquid 4 He is accurately known from 
previous measurements on MARI using a high resolution 
setting. 26 Using this J{Q, y) as input, the present broader 
instrument resolution can be determined directly by re- 
producing the observed resolution broadened J(Q,y) of 
pure He {x = 0). The resolution determined in this way 
at Q = 27 A -1 is shown in Fig. [T] 



Standard procedures were employed to convert the 
raw neutron scattering data from TOF and scattered 
intensity to energy transfer and the dynamic structure 
factor, S(Q,w). A summary of these conversions has 
been given by Andersen et al.^L The data was then con- 
verted to the y scaling energy transfer variable and to 
J{QiV) = vr.S(Q,u>). The data were analyzed in detail 
for wave vectors 26 < Q < 29 A -1 . 

The background arising from the empty cell was mea- 
sured separately and found to be a smooth function of 
scattering angle and TOF with no microstructure. At a 
given Q value, the measured background was well repre- 
sented by a smooth quadratic function in y, ay 2 + by + c, 
where a, b and c were determined by a least square fit 
to the measured background. An example of this back- 
ground is shown in Fig. |2]as a dotted line. This measured 
background function was subtracted from the scattering 
from the helium plus cell to obtain the net measured 
scattering intensity from liquid helium. 

The analysis consists of representing the net J(Q, y) as 
a sum of models for the scattering 4 He and 3 He. A linear 
background of variable height and slope was also included 
to allow for any error in the background determination 
or shielding of the background by the sample. This lin- 
ear background was always very small. Additional error 
in model parameters for 4 He and 3 He arising from incor- 
porating the linear background is included in the quoted 
parameter errors below. The data was y scaled so that 
the 4 He peak is centered at y = using the 4 He mass. 
This means that the 3 He peak component of the model 
function must be multiplied by M — m^jm^ to regain 
the correct magnitude. An example of the total scat- 
tered intensity y scaled and of models 4 He and 3 He plus 
the measured background is shown in Fig. [21 
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FIG. 2: Total observed scattering intensity (solid circles) showing 
J(Q> y) f° r a 3 He- 4 He mixture at 10% 3 He concentration and T = 
2.5 K plus background. Dotted line is a separate measurement of 
the background. The data is presented with the 4 He peak centered 
at j/4 = 0. The error bars are the size of the circles. 



III. DATA ANALYSIS 

Our goal is to determine the condensate fraction and 
to get as much information as possible about the He 
and 3 He momentum distributions in the mixture. The 
data is good enough to determine at most two free fit- 
ting parameters for each peak. To describe the 4 He peak, 
we follow the method used for the pure liquid 4 He. In 
superfluid 4 He, we use the convolution approach (CA) 
in which a FS broadening function R(Q, y) is convoluted 
with a model Jja(h) as in (JSJ. In normal 4 He, we use the 
additive approach (AA), described below, and the CA. 
The condensate fraction and the width (~ kinetic energy) 
of the peak are the free parameters. The FS function and 
the shape of n(k) are assumed to be the same as in pure 
4 He. The comparison between normal and superfluid liq- 
uid 4 He in Fig. 2] shows that the peak height is higher 
at y = for the superfluid phase due the presence of 
a condensate and a left-right asymmetry introduced by 
the condensate term noR(Q,y) in J(Q,y). To describe 
the 3 He peak, we construct a model of a Fermi momen- 
tum distribution n(k) and compare the corresponding 
resolution broadened Jia (y) in (0) with the data (no FS 
effects). The parameters in the 3 He n(k) are the step 
height at the Fermi surface and the length of the high 
momentum tail. 

There is a small overlap of the net 3 He and 4 He peaks. 
Fig. [21 shows that the degree of overlap is very small 
indeed. The area of overlap of the models is typically 
1-2% of the total area under the two models. We found 
that the parameters in the model fit to one peak were 
independent of the model used for the second peak within 
the error quoted. This was tested by using various models 
to represent the second peak. 



A. 4 He response 

1. Convolution Approach 

The convolution approach (CA) is most useful when 
the momentum distribution n(k), and thus JiA{y) in J2J, 
are narrow in y relative to R{Q, y). Equally, it is useful if 
n(k) contains a component, such as a condensate, that is 
narrow relative R(Q, y). In this case, the observed width 
of the narrow component in J(Q,y) is set by R(Q,y). 
In pure liquid 4 He, the FWHM of R(Q, y) is approxi- 
mately 1 A -1 . This is broad compared to a condensate 
component but narrower than the remainder of n(k). 

To determine the condensate fraction, it is convenient 
to separate the state that is macroscopically occupied 
from the regular uncondcnscd states. In a uniform liquid, 
the natural orbitals are plane wave, momentum states. 
The orbital containing the condensate is the k = state. 
This state contributes a term nn<5(k) to n(k) where no = 
Nq/N is the fraction of particles in the condensate. The 
regular, uncondensed states are the k ^ states and 
we denote their contribution to n(k) as n*(k). In an 
interacting Bose liquid, Bosons can scatter into and out 
of the condensate from the k ^ states. This leads to 
a coupling between the condensate and the k ^ states 
and a terr^Z^ 
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in n(k). The coupling is strongest for the low k states 
and the expression in the square bracket of (@J is derived 
for and valid for the low k limit. We have multiplied it 
by a Gaussian to cut off no/ (k) at higher k as done for 
pure2& 4 He with k c fixed at k c = 0.5 A -1 . This term 
is highly localized around k = and when Final State 
broadening is included no/(k) cannot be distinguished 
experimentally from no<5(k). The model n(k) we use for 
the 4 He component therefore has three terms2i, 



n(k) = no[<J(k) + /(k)]+Am*(k) 



(5) 



Ai is a constant chosen by normalization J <ik?i(k) = 1. 

To implement this model, we Fourier transform n(k) 
to obtain the corresponding One-Body-Density-Matrix 
(OBDM), n(r)=f dke lk r n(k). J IA (s) is the OBDM for 
displacements r = sQ parallel to the scattering wavevec- 
tor Q. Jia(s) corresponding to (J3J) is obtained by Fourier 
transforming JSJ to obtain n(r) and using Jia(s) = n(s) 
giving, 



Jia(s) = n(s) = n [l + f(s)] + Am*{s) 



(6) 



The term [1 + f(s)] is long range in s. The n*(s) is 
short range in s and we represent it by2& 



n*(s) = exp 
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(7) 
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where 0,%, 0,4 and ae are parameters (cumulants). This 
representation is useful when n*(s) is at least approxi- 
mately a Gaussian, plus small corrections. In going from 
iJSJ to |JBJ, we actually need to transform /(k) only, and 
once only since /(k) has no free parameters. Given /(s), 
our model may be viewed as a model for Jja{s) in 

The model Jia(s) is multiplied by the FS function 
R{Q, s) and the product is Fourier transformed to obtain 
J(Q, y) as in 10. In the present fit to the 4 He peak, we 
used the pure 4 He FS function at all 3 He concentrations, 



R{Q, s) = exp 



0: 



3 3 ^5 5 P6 6 

3l S "IT ~ 6! S 



As 



(8) 



with parameters (3% = 0,3/ XQ, (3$ = a^/iXQ) 3 and 
= &64/(AQ) 2 with (13,052 and 0,54 set at their pure 
4 He values^ We also set 0:4 and 6tQ in n*(s) at their 
pure 4 He values^ The J(Q,y) therefore has only two 
free parameters, ai and uq. 

The above procedures for determining the shape of 
the momentum distribution of pure liquid 4 He, and ex- 
tracting the condensate fraction are well established, 
and have been the subject of intense activity over the 
years, 21 . 26 ! 47 . 52 



2. Additive Approach 

When the momentum distribution is broad, FS effects 
are relatively less important. In this case, the exponential 
in R(Q, s) of JSJl can be expanded and FS effects retained 
as additive corrections to the IA. Similarly, if the devia- 
tions of n*(s) from a Gaussian are not too large, n*(s) 
in can be expanded and deviations from a Gaussian 
n*(s) retained as additive corrections. This expansion 
of n*(s) is the same as the Gauss- Hermite expansion in- 
troduced by Sears* 2 ^ These expansions lead to the Addi- 
tive Approach (AA), in which J(Q,s) is represented as 
a Gaussian IA, Jja{s), plus corrections for deviations of 
n*(s) from a Gaussian and for FS effects &$M££°, 



J(Q, s) = J IA (s) + Ji(Q, s) + J 2 {Q, s) 



(9) 



where 



Jia(s) = (27ra 2 ) 



-1/2 e ~s 2 /2a 2 



M s ) = y^3S 3 JlA(s) 
J2(s) = ^4S 4 J/a(s) 



(10) 
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FIG. 3: Schematic diagram of the Fermi model momentum dis- 
tribution used in this analysis. The parameters are defined in the 
text. 



B. 3 He Response 

The 3 He momentum distribution is not well approx- 
imated by a Gaussian. Thus we do not expect the fit 
of a Gaussian or the A A to the 3 He peak to reveal the 
physics well. To proceed, we constructed a simple model 
of the 3 He n(k) consisting of a step of height Z at the 
Fermi surface plus an exponential high momentum tail. 
We calculate the Jia{v) for this n(k) using J2J and fit 
the convoluted Jia(u) directly to the data. 

Specifically, the model n(k) is, 



n(k) 



H 

V F 



-r- < 1 
r- > 1 



(11) 



where kp is the Fermi momentum and Vp — ^§-kp the 
volume of the Fermi sphere. Assuming that the volume 
occupied by 3 He in the mixtures exceeds that occupied 
by 4 He by about 28%, as suggested by the dielectric 
measurements of Edwards et al^, the Uf can be ap- 
proximated using the expression kp = (37r 2 px) 1 ^ where 
p is the average density of the mixture at SVP and x 
the 3 He concentration. The kp values are listed in Ta- 
ble [I] The model has three parameters, H, Hp and 
77 where 77 determines the length of the high momentum 
tail and Z = H - Hp (See Fig. 0). Normalizat ion of 



and P3 = CI3/XQ and pi — 0:4 + a^KXQ) 2 . The 6:4 is the 
leading deviation from a Gaussian and 03 and 0,4 are FS 
terms. In fitting this function to the 4 He peak in normal 
mixtures, we fixed 03 and a.4 at their pure 4 He values 
(S4 = 0). The J(Q, y) then has only one free parameter, 
CI- 



TABLE I: Fermi momentum kp(x) = (37r 2 px) 1 / 3 as a function of 
the 3 He concentration x where p is the average mixture density at 
SVP. 
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n(k), J dkn(k) = 1, gives 

3H F ( 2 2 \ 

H+ + _ + l (12 

rj \r) 2 7] J 

which can be used to eliminate one parameter, say H. 
Similar models of n(k) are discussed by Carlson et a/4* 
and used by Azuah et al£& and Mazzanti et aim* for pure 
liquid 3 He. The IA given by corresponding to (|llfl is, 

7 , , fw(*(w 2 ) + ^(i+±)) y'< 1 

(13) 

where y' = y/kp. 

If the liquid temperature T is near or above the Fermi 
temperature, Tf = 2Tr f» ks kp , there may be some thermal 
broadening of the Fermi liquid n(k). 20 This is certainly 
the case for a Fermi gas for which n(k) is a step func- 
tion (Z = 1) at T = K. However, in the liquid model 
above, interaction has already reduced the magnitude of 
the step at kp (Z < 1). Thus thermal broadening may be 
relatively less important in a strongly interacting liquid. 
This is discussed more fully in section llVl 

To estimate Tf, we note that the effective mass m* 
depends strongly on the pressure but little on the 3 He 
concentration^ We take m* — 2.3m as found by theoret- 
ical calculationsi2^£, independent of concentration. In 
the range x = 5-20 %, the corresponding Tp is 0.3-0.8 K. 
This is significantly less than the normal liquid temper- 
ature T = 2.5 K and comparable to the superfluid tem- 
perature T — 0.4 K used here. Our n(k) extracted from 
data will represent n(k) at these temperatures. However, 
we did not observe any difference between n(k) at T = 
0.4 K and 2.5 K. This suggests that the rounding of the 
small step Z at kp by FS effects and the instrument res- 
olution function is comparable or greater than thermal 
broadening. 



IV. RESULTS 

In this section, we present our neutron scattering data 
and the results for the 4 He atom kinetic energy, the 3 He 
momentum distribution and the 4 He condensate fraction 
obtained from analyzing the data. Since the absorption 
cross section of 3 He is so large, the net scattered neutron 
intensity from 3 He- 4 He mixtures is weak. For this reason, 
a broad energy resolution setting was used on the MARI 
spectrometer, as shown in Fig. ^ Also, the statistical 
precision of the data is not high. As a result, we were 
able to determine uniquely only one or two parameters 
in model fits to the data. Specifically, we found that 
the parameter no that provides the condensate fraction 
and the parameter cxi that sets the kinetic energy, K4 = 
3h 2 ct2/2m, were correlated. For this reason, we begin 
with normal 4 He where uq = 0. 



A. Normal phase 

1. 4 He kinetic energy 

The kinetic energy K4 of 4 He in the normal liquid at 
T = 2.5 K is shown in Fig. 0|as a function of 3 He concen- 
tration. K4 is determined from the 4 He peak using two 
fitting procedures, the convolution approach (CA) and 
the additive approach (AA). In the CA, we use the Final 
State function that was determined^ previously for pure 
4 He unchanged at all 3 He concentrations. The only free 
parameter in the present fit was 612 ■ Similarly, in the AA 
the FS parameters were set at their pure 4 He values 2 ^ 
and only 0,2 was free. The parameters held fixed in the 
CA are summarized in Table ITT1 

Fig. 21 shows that the CA and AA procedures give 
identical kinetic energies for pure 4 He (x = 0). This 
value agrees with our previous^ 4 He determination 
K4 =16.3±0.3 K and with the value found by Sen- 
esi et al.m^ As the 3 He concentration is increased, the 
two methods continue to agree within error but the CA 
method gives a marginally lower K4. The K4 decreases 
only very little with increasing x. Azuah et alm& and 
Senesi et alM find that K4 decreases more sharply with 
increasing x. A decrease is anticipated since the density 
decreases as x increases and in pure fluids K4 decreases 
as density decreases. When we fit a Gaussian to the 4 He 
peak, we find a K4 that is independent of x. We return to 
this point in section TlV CI and compare with calculations 
in section IVI 

The K4 values in Fig. 0]are averages obtained from fits 
to data at specific Q values. The variation of K4 with Q 
is shown in Fig. [S] The aim is to display the statistical 
precision of the data. 



2. 3 He momentum distribution 

In this section, our goal is to learn as much as possible 
about the 3 He momentum distribution. We do this by fit- 
ting the JiA(y) (|13|) obtained by substituting the model 
n(k) given by l|ll|) and shown in Fig. |3J into 10) • The 
model n(k) has two parameters, Hp giving the height 
of the high momentum tail at kp and 77 which sets the 
length of the tail. H is determined by normalization. 



TABLE II: Fitting variables used in the CA. In the normal phase, 
only Q2 is determined from a fit to data and in the superfluid phase 
only no is determined. The values of fc c , j3 n and a„ for n > 3 were 
all kept fixed at their values in the pure 4 He case, as evaluated by 
the precise measurements of Clyde et al.ii 
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FIG. 4: 4 He kinetic energy [Jf 4 ] at T = 2.5 K as a function of 3 He 
concentration, x. The closed squares are the 4 He kinetic energies 
obtained from the Additive Approach (AA) and the open circles 
obtained from the Convolution Approach (CA)™i 



The data in 3 He peak region was precise enough to de- 
termine one parameter well with an estimate of a second 
parameter. Since the tail of n(k) is of specific interest, we 
also set the tail parameter at specific values e.g. no tail 
(H = 1), tail that corresponds to the calculated n(k) 39 
(77 = 0.8) and fitted for Hp to see how good a fit could 
be obtained. We subsequently obtained best fit values 
of 77. Once, 77 is determined, the model kinetic energy is 
obtained as K3 = J c?k/c 2 n(k) giving, 



3h 2 k 2 F 

2m,3 5 



H 



5H F 4! 

V 4 
?7 77* 



4! 



12 

1^2 



4 
V 



r 



(14) 

where ^— = 8.08 K A 2 . For typical values of 77, more 
than 95% of K3 arises from the tail in n(k). Also, in 
(|14fl . K-j, is very sensitive to 77. Thus, unless 77 can be 
accurately determined, K3 is not a good single parameter 
to characterize the data. 

Fig. shows fits to the 3 He peak with no high mo- 
mentum tail (Z = H = I, Hp = 0) and with the tail 
parameters set at 77 = 0.8, the value obtained in a fit to 
the high k portion of the tail calculated by Boronat et 
al. Clearly a Fermi step function with no tail cannot 
reproduce the data. In contrast, an n(k) with a high en- 
ergy tail that reproduces the calculated tail (see Fig. EjJ 
fits the data. In this sense, the calculated n(k) and the 
data are entirely consistent. There is no disagreement 
between theory and experiment. 

We determined an observed value of 77. The best fit 
value of 77 obtained varied somewhat with the data set 
(e.g Q value) considered reflecting the statistical pre- 
cision of the data. Fig. shows fits for 77 = 0.76 
and 0.94 which provide equally good fits at different 
Q's. The data is precise enough to determine that 77 
lies in the range 0.7 < 77 < 1.0, i.e 77 = 0.85±0.15 and 
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FIG. 5: Top: Additive Approach (AA) fit to data in the normal 
phase (T = 2.5 K) at several 3 He concentrations showing the Q 
dependence of K4. Bottom: Convolution Approach (CA) fit to 
data in the normal phase. 



Z = H — Hf = 0.05 ± o of . This correspond to a very 
strongly interacting Fermi liquid. Fig. [S] shows a com- 
parison between our observed n(k) for a 5% mixture at 
T = 2.5 K and a calculated^ n(k) at 6.6%. 



B. Superfluid phase 

As noted, we found that the parameter cb and the 
condensate fraction no were correlated. For example, in- 
creasing no and decreasing both lead to a narrower 
peak that cannot be distinguished within the present pre- 
cision of the data. In pure 4 He, a.^ was found to be the 
same in normal (T = 2.5 K) and superfluid (T = 0.5 
K) 4 He. Thus in the mixtures we set ct2 in the super- 
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FIG. 6: Fits to the 3 He peak (solid lines) in the normal helium mixture data at x = 5% and Q = 29 A -1 . Data are open 
circles with error bars. The error bars are the size of the circles except where shown. The left side shows a fit assuming 
the 3 He rt(k) is a Fermi step function (i.e n(k) has no high momentum tail (H =1 )) a nd the right side shows a fit to data 
with the model n(k) shown in Fig. |J]that reproduces the calculated tail of Ref. I.'>!l . The fit with no tail (left) is poor 
while the fit incorporating the calculated tail (right) is good. 




FIG. 7: Fits to the 3 He peak (solid lines) in the normal helium mixture data (open circles) at x =5% and Q = 28.5 
A" 1 . The error bars are the size of the circles except where shown. The graphs show that it is possible to get good fits 
to the same data with two different sets o f pa rameters. This illustrates the precision in which the n parameter can be 
determined. The calculated tail from Ref. l3fll corresponds to r\ = 0.8. 



fluid phase at the normal phase value determined above 
at each x. 

1. Condensate fraction 

To determine no (x) , we first perform fits of the C A to 
data for each Q at a given concentration x. The average 
value over all Q gives the corresponding no (x) . The varia- 
tion of the condensate fraction with Q is shown in Fig. [5] 
The variation reflects the statistical precision of the data. 
Fig. 1101 shows the resulting x dependence. The helium 
mixtures data at x = 20% was taken at a slightly higher 
temperature (T = 1.3 K) than the lower concentrations 



data. To be able to make a consistent comparison of the 
results on n Q , we made a temperature correction^ to n 
at x = 20% using n (T) = n [l- (^) 7 ] where T\ - 2.17 
and 7 = 5.5 to get its equivalent value at T = 0.4 K. The 
corrected n Q is shown in Fig. ^|along with the results at 
lower 3 He concentrations. We find that no (a;) increases 
slightly above the pure 4 He value with increasing x. Fig. 
1101 represents our final results for no which are discussed 
in section fvl 
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FIG. 8: Comparison of the observed 3 He atom momentum distri- 
bution n(k) (solid line) at x =5% with the calculated n(k) from 
Ref. 15s! at x = 6.6% (dotted line). 



2. 



i He kinetic energy 



The 4 He kinetic energy in the superfluid phase can be 
obtained using the form of the momentum distribution 
in ijjjj and assuming that the parameter a.i is the same in 
both normal and superfluid phases. Normalizing the He 
model momentum distribution, as discussed above, leads 
to n [l + If] + At = 1 with I f ~ 0.25. This result can 
be obtained from © as n(s — 0) where If = f(s = 0) 
and n*(s = 0) = 1. Since n Q is a function of x, the 
normalizing constant Ai is also a function of x, Ai(x) = 
1 — 1.25tt,o(.t). The kinetic energy K4 is proportional 
to the second moment of n(k). Only n*(k) contributes 
significantly to the second moment so that from JSJ and 
with a2 the same in the normal and superfluid phases we 
have K^(S) = A%(x)K4(N). In this model, the kinetic 
energy in the superfluid phase (S) is reduced below that 
in the normal phase (N) entirely by BEC. We found the 
kinetic energy to be about 10% lower in the superfluid 
phase than in the normal phase increasing somewhat with 
x as shown in Fig. 1111 



3. 3 He momentum distribution 

Fig. 1121 shows two fits to the 3 He peak in the super- 
fluid mixture at T = 0.4 K. In one fit, the tail parameter 
of the model n(k) is set at 77 = 0.8, the value that re- 
produces the calculated tail. The same model (77 = 0.8) 
was compared with data for T = 2.5 K in Fig. [fj] The 
agreement with the data for 77 = 0.8 at T = 0.4 K and 
2.5 K is indistinguishable. This shows firstly that there is 
no observable temperature dependence of the 3 He peak, 
between T = 0.4 K and T — 2.5 K which corresponds to 
T ~ Tp and T ~ 6Tp, respectively. 
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FIG. 9: The 4 He condensate fraction (closed squares), no, as a 
function of momentum transfer Q. The error bars represent statis- 
tical errors obtained from fits to data. The solid lines are guides 
to the eye. 



Somewhat surprisingly, there is no observable temper- 
ature dependence of the 3 He momentum distribution al- 
though the two temperatures are significantly different. 
We return to this point in the discussion section. Sec- 
ondly, there is again no disagreement between theory (for 
T = 0K) and experiment. 

As at T = 2.5 K, a range of r\ values provided best fits 
to data depending upon the Q value considered. Fig. IT21 
shows an example at Q = 27.5 A -1 where the best fit is 
obtained for r\ = 1.0. At T — 0.4 K, we find r\ lies in the 
range 0.7< r] < 1.0 or 77 = 0.85±0.15 as found at T = 2.5 
K. 



C. Alternate methods of analysis 

In this section, we test the sensitivity of the results to 
the method of analysis by analyzing the data in different 
ways. 
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FIG. 10: The 4 He condensate fraction (closed squares), no, as a 
function of 3 He concentration. The error bars are standard devia- 
tions obtained from a linear least-squares fit to the Q dependence 
of no- The condensate fraction increases with increasing x from 
7.25± 0.75% at x = to 11±3 at x = 15-20%. The dashed line is 
a guide to the eye. 
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FIG. 11: Kinetic energy of 4 He atoms, K4, as a function of 
concentration. The superfluid K4 were obtained from the 
normal K4 assuming the condensate does not contribute to 
the kinetic energy. 



1. Kinetic energy assuming a Gaussian 

A straightforward method of analysis is to fit Gaussian 
functions to both the 3 He and 4 He peaks as employed by 
Azuah et al.w& This method, when applied to our data 
yields a 3 He kinetic energies that agree with those previ- 
ously obtaine d 29 i 31 where a Gaussian or a Gaussian-like 
fit was made. Senesi et alM^L also included the leading 
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FIG. 12: Fits to the superfluid helium mixture data at a 3 He 
concentration x = 5%. The 4 He peak is fitted using the CA method 
and the 3 He peak assuming a Fermi n(k) having a variable tail 
length n. The open circles are data points and the solid lines are 
fits to data. The error bars are the size of the circles except where 
shown. 



term beyond a Gaussian in a Gauss-Hermite expansion 
of J(y) to analyze their data. However, they found that 
the inclusion of higher order terms did not result in sig- 
nificant improvements in the fits. Both Azuah et al. and 
Senesi et al. found K4 values that decrease with increas- 
ing 3 He concentration. 

In Fig. El we show the K4 values obtained by fitting a 
Gaussian J(Q, y) to our data compared with the reported 
values by Azuah et alm& and by Senesi et al.m^ Firstly, 
for pure 4 He (x = 0) we find a K4 that lies approximately 
10% below the CA values. A smaller K4 is expected for 
a Gaussian. In the CA, FS effects cut off the wings of 
the peak and allow a somewhat broader ?i(k) to fit the 
data. Senesi et al. took data at Q ~ 100 A -1 where 
FS effects are probably negligible. Thus we expect their 
values to agree with our CA values. Most importantly, 
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FIG. 13: Observed K4 (open circles) assuming a Gaussian fit to 
the 4 He J(Q,y). The measured values by Azuah et al. (open 
triangles) 29 and Senesi et al. (open squares)^™ are shown for com- 
parison. The dashed lines are guides to the eye. 



FIG. 14: Condensate fraction, no, obtained from present data 
assuming the kinetic energies measured by Azuah et al. and Senesi 
et al£L. The no decreases with increasing x. This is illustrates the 
correlation between no and K4. 



using a Gaussian fit we find a K4 that is independent of 
x. Thus while our CA, the Azuah et al. and the Senesi et 
al. K4 agree for pure 4 He value (x = 0), the data itself 
appears to disagree on the x dependence of K 4 . Most 
theoretical calculations also find that K4 decreases with 
increasing x as shown in Fig. 115b . 

2. Condensate assuming concentration dependent kinetic 
energies 

Since there is some correlation between the c\i param- 
eter (if 4) and no, we determined no using the values of 
CY2 corresponding to the measured 4 He kinetic energies 
by Senesi et al.~ The resulting values of no are shown 
in Fig. 1141 We find no decreases with increasing x, when 
the x dependent K4 values of Senesi et al. and Azuah et 
al. are used. This result is expected since a narrower 4 He 
peak is obtained by either increasing no or decreasing 0.2- 
If c<2 (~ K4) decreases, a smaller no is needed to obtain 
a good fit. Thus a K4 that decreases significantly with x 
implies an no that also decreases with x in disagreement 
with theory. 

V. DISCUSSION AND CONCLUSION 

A. 4 He response 

In pure normal 4 He, we find a kinetic energy [K4 = 
16.1 ± 0.3 K) which agrees well with previous neu- 
tron scattering data™ (K4 = 16.3 ± 0.3 K ) and with 
calculations^ (K4 = 15.41 K). In normal helium mix- 
tures, we find a K4 that decreases somewhat with in- 



creasing 3 He concentration. This suggests that the lo- 
cal environment of individual 4 He atoms remains largely 
unchanged^ with the addition of 3 He atoms at low x. 
This finding represents the chief difference between our 
results and early theoretical treatments" and previ- 
ous experiments 29,31 which find that K4 decreases signif- 
icantly with increasing x. Theoretical and experimental 
values of K4 are compared in Fig. 115b . The present K4 
and a more recent Path Integral Monte Carlo (PIMC) 
calculation^ of K4 which show a substantial agreement 
on the x dependence are compared separately in Fig. 115b . 
Clearly, there remain some theoretical and experimental 
differences to be resolved. The present K4 in the super- 
fluid phase was obtained from the normal phase value 
assuming that the condensate does not contribute to the 
kinetic energy. 

To determine the condensate fraction, no, in helium 
mixtures, we introduced a model n(k) which has a con- 
densate term and fitted the model to the neutron scat- 
tering data. In the model, the second moment of n(k) for 
the finite k states was assumed to be the same as in the 
normal phase, as we found previously 2 ^ for pure liquid 
4 He. The higher moments and the Final State function 
were assumed to be the same as in pure superfluid he- 
lium. 

In mixtures we find a small increase in no with in- 
creasing 3 He concentration, to no = 11 ± 3% at x = 
15-20 %. Fig. 1161 shows our observed no as a function 
of x compared theoretical values. Our results agree with 
the theoretical values within experimental error. Wang 
and Sokol 34 measured no at one concentration reporting 
a value of no = 18% at x = 10%. This is higher than 
the present measured and calculated values. The con- 
centration dependence of no may be roughly estimated 
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FIG. 15: (a): 4 He kinetic energy (K4) as a function of 3 He concentration, theory and experiment. Calculations are VMG™ (solid 
circles), HNC using a Jastrow function™, (triangles), HNC including triplet correlations^? (solid squares). Experiments are Azuah et alm& 
(up-pointing triangles), Wang and Soko 134 (down-pointing triangle), Senesi et aldil (open squares) and present at T = 2.5 K (open circles). 
Dashed lines are fits to experiments, (b): Present observed K4 at T = 2.5 K (open circles) as a function of 3 He concentration compared 
with a PIMC calculation™ of K4 at T = 2.0 K (solid stars). Dashed lines are guides to eye. 



assuming that the effect of the 3 He is simply to change 
the volume available to the 4 He (average density approx- 
imation) and using the density dependence of no in pure 
4 HeS Essentially, the molar volume of the mixture at 
low concentrations is approximately^ v ~ V4 (1 + ax) 
where V4 is the 4 He molar volume at SVP and a the ex- 
cess of volume occupied by 3 He (~ 0.28). This gives a 
mixture density which decreases with x. On the other 
hand, the condensate fraction, no, in pure liquid He 
increases22i££ with decreasing density. This yields an no 
which increases very slightly with x, from no ~ 7% (x — 
0) to n ~ 9% (x = 20%). This simple calculation sug- 
gests that no should increase little with concentration, in 
qualitative agreement with our experiment. 
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We find that a model 3 He momentum distribution n(k) 
with a small step, Z, at the Fermi surface (k — kp) and 
a substantial tail at higher k characteristic of a strongly 
interacting Fermi liquid reproduces the observed 3 He re- 
coil peak well. For a gas of non-interacting Fermions at 
T = K, n(k) is a simple step function (Z —1); n(k) = 1 
for k < kp and n(k) = for k > kp. This case is repro- 
duced by our model n(k) in Eq. (|ll(l for parameters H = 
1, Hp = (Z = H - Hp = 1, no tail in n(k)). When 
Fermions interact, some Fermions move from states be- 
low kp to states above kp to minimize the total energy, 
potential plus kinetic energy. The step at kp is reduced 
Z < 1 and n(k) develops a tail at higher k. This is repro- 
duced in Eq. for parameters H — Hp = Z < 1 and 
77 > 0. The stronger the interaction, the smaller is Z and 



FIG. 16: The condensate fraction (no) as a function of 3 He con- 
centration, theory and experiment. Calculations are the open sym- 
bols; VMG™ (circles), HNC (Lennard- Jones )Sm (down-pointing tri- 
angles), HNC (triplet correlations )22 (squares), HNC (Jastrow)ffi 
(up-pointing triangles), HNC (Aziz potential) 3 ?, (right-pointing tri- 
angles). Measurements are Wang and Soko! 34 (solid circle) and 
present (solid squares). 



the longer is the tail (smaller 77). In pure liquid 3 He a 
recent diffusion Monte Carlo (DMC) calculation 4 ^ finds 
Z = 0.236 and a tail in n(k) that is well represented by 
an exponential (times algebraic terms) in k. This n(k) 
agrees well with earlier DMC results 

Pure liquid 3 He is regarded as a strongly interacting 
Fermi liquid, more strongly interacting than nuclear mat- 
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ter, and Z ~ 0.25 represents a large reduction of the step 
height at kp arising from interaction. As noted above, 
Boronat et al£^ calculate an n(k) obtaining Z = 0.09 and 
a tail that is also well fitted with an exponential with 77 = 
0.8 for a 3 He- 4 He mixture of x — 6.6 % 3 He. Apparently, 
3 He in a dilute mixture is more strongly interacting than 
in pure 3 He. Our observed, finite temperature n(k) has 
Z = 0.05 ± Ej^ 4 and r] = 0.85 ± 0.15 at x = 5 %. It agrees 
with the calculated n{k) for T = K. Z = 0.05 ± g;Jf 
represents a very small step at kp . 

We have not included Final State (FS) effects in our 
analysis of the 3 He recoil peak. The impact of FS ef- 
fects depends critically on the width and shape of the 
intrinsic Jia{d) arising from n(k). If Jia{v) is narrow 
relative to the FS broadening function R(Q, y) or has a 
narrow component or sharp feature, then the FS effects 
serve to broaden the narrow component to the width of 
the FS function. The classic example is the condensate 
component in superfluid 4 He. However, if Jja(]j) is broad 
relative to the FS function, then the FS function has little 
further broadening impact. For example, the second mo- 
ment of the FS function is zero. In this event the FS ef- 
fects can be well represented by a moment expansion 4 ^^ 
in which the leading term is the third moment. This term 
serves to introduce an asymmetry into J(Q,y) but will 
not modify its width or whether J(Q, y) has significant 
tails at larger y arising from tails in n(k). 

The FWHM of R(Q, y) at Q w 25 A" 1 is in the range 
0.7-1.0 A -1 in pure liquid 3 He, in 3 He- 4 He mixtures and 
in pure liquid 4 He. 26 i 44 i 49 The FWHM of Ji A (y) of 3 He 
in 3 He- 4 He mixtures is approximately 2.0 A -1 , largest 
in dilute 3 He mixtures. We expect FS effects firstly to 
spread out the discontinuity in Jia(jj) arising from the 
step of height Z in n(k). It will be spread out over a 
width y ~ 0.7 A -1 (i.e. of order of 2kp). The present 
instrument resolution function, shown in Fig. ^ which is 
broader than the FS function at Q ~ 25 A -1 , will further 
spread out the discontinuity. Parenthetically, it is prob- 
ably largely because of the resolution and FS broadening 
that we do not observe any apparent broadening of n(k) 
arising from temperature. Otherwise, the FS function is 
narrow relative to the width of Jia(d) arising from n(k). 
Thus we do not expect FS effects to change the overall 
width of J(Q,y) (i.e. the apparent width of n(k) and 
whether it has tails or not). However, including FS ef- 
fects should improve the shape of J(Q, y) as it has been 
done^ in pure liquid 3 He. Thus we do not expect FS 
effect to alter our basic conclusions about n(k). 

We note that one cannot use the FS function for pure 
3 He for mixtures. Essentially, in dilute 3 He mixtures the 
Fermi statistics for only the 3 He atoms should be in- 
cluded rather than for all atoms as in pure 3 He. We plan 
to incorporate FS effects in a future publication. 

We have strictly used a model n(k) valid for T = K 
since the step Z at kp is not rounded by thermal effects 
to fit data taken at T = 0.4 K ~ Tp and T = 2.5 K 
~ 6Tp. Stated differently, we obtain an n(k) at two finite 
temperatures from the data. These n(k) are expressed 



in the form of a T = K n(k). The n(k) at the two 
temperatures are the same within experimental error. 

For a Fermi gas, the thermal broadening of n(k) is 
well known. At T ~ Tp the broadening is significant 
with n(k = 0) = 0.73, n(k F ) = 0.5, n(2k F ) = 0.23. It 
could be said that there is a rounded step of Z ~ 0.5 in 
ri(k) over a width of k ~ 2kp and n(k) develops a short 
tail. This thermal broadening is less than the resolution 
and FS broadening. The thermal broadening of n(k) for 
an interacting Fermi liquid is not known. However, in a 
strongly interacting Fermi liquid Z is already small and 
n(k) already has a long tail. Hence thermal broadening 
is expected to be relatively less important in an interact- 
ing liquid. For example, a higher temperature might be 
required before the thermal rounding (e.g. effective ther- 
mal reduction of Z) becomes significant relative to the 
'interaction' reduction. Apparently, this rounding is not 
significant or critical compared to the rounding arising 
from FS effects or the instrumental broadening within 
current error. The model parameters Z = 0.05 ±0 01 and 
77 = 0.85 ± 0.15 obtained are consistent with theory for 
T = K. Stated differently, the calculated n(k) for T = 
K reproduces our experimental data. Our observed Z 
may be somewhat smaller than the calculated value be- 
cause of resolution and thermal effects. The theory and 
experiment are therefore consistent in finding that 3 He in 
dilute 3 He- 4 He mixtures is a strongly interacting Fermi 
liquid. This strong interaction must arise at least in part 
from interaction with the 4 He component of the liquid in 
dilute mixtures. 

Mazzanti and co-workers 44 have compared their calcu- 
lated J(Q,y) for pure liquid 3 He with that observed by 
Azuah et al. at T = 1.4 K (T ~ Tp) and Q = 19.4 A" 1 
in much the same way we have done here. The input to 
their calculated J(Q,y) is n(k) for T = K calculated 
using DMC. This yields Z = 0.236 as noted above. They 
represented their calculated n(k) by an analytic function 
that had an exponential tail times algebraic terms. The 
second input is a Final State broadening function valid 
for a pure Fermi liquid. The Final State broadening func- 
tion R(Q, y) has a FWHM in y of approximately 0.7 A -1 
compared with the total FWHM of J(Q, y) of 2.0 A" 1 so 
that R{Q, y) has only a modest impact on J(Q, y). They 
find their J(Q,y) reproduces the observed J(Q,y) well 
and conclude that here is agreement between theory and 
experiment. In a similar way we have used a model n(k) 
with parameters (rather than a calculated n(k)) repre- 
sented by an analytic function with an exponential tail. 
We find a J{Q, y) (without FS effects) which reproduces 
experiment within statistical error. The model n(k) is 
consistent with n(k) calculated for mixtures^ We con- 
clude in a similar way that there is no disagreement be- 
tween theory and experiment for mixtures as well as in 
pure 3 He. 

In summary, we have measured the momentum distri- 
bution of 3 He- 4 He mixtures for 3 He concentrations x be- 
tween and 20%. We find uq increases from 7.25±0.75% 
in pure 4 He (x = 0%) to 11 ±3% for mixtures with 
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x = 15 — 20%, in agreement with theoretical calculations. 
The 4 He kinetic energy is found to decrease slightly with 
3 He concentration. The model 3 He n(k) that reproduces 
our data at x = 5% agrees with the calculated n(k) at 
x = 6.6% removing a previous apparent disagreement 
between theory and experiment based on comparing 3 He 
kinetic energies. 
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We present neutron scattering measurements of the momentum distribution of liquid 3 He- 4 He 
mixtures. The experiments were performed at wavevectors Q, 26 < Q < 29 A -1 , on the MARI time- 
of-fiight spectrometer at the ISIS Facility, Rutherford Appleton Laboratory, a spallation neutron 
source. Mixtures with 3 He concentrations x between and 20% were investigated both in the 
superfluid and normal phases. From the data, we extract the Bose-Einstein condensate fraction no 
and the momentum distributions of 3 He and 4 He atoms. We find that no increases somewhat above 
the pure 4 He value when 3 He is added; e.g from no = (7.25 ± 0.75)% at x — to (11 ± 3)% at x = 
15-20%. This agrees with predictions but is less than the only previous measurement. We find a 4 He 
kinetic energy K4 for pure 4 He that agrees with previous determinations. K4 decreases somewhat 
with increasing 3 He concentration, less than observed previously and found in early calculations but 
in agreement with a more recent Monte Carlo calculation. The 3 He response is not well reproduced 
by a Fermi gas momentum distribution, n(k). Rather an n(k) having a small step height at the 
Fermi surface and a substantial high momentum tail characteristic of a strongly interacting Fermi 
liquid provides a good fit. This n(k) is consistent with calculated n(k). Thus agreement between 
theory and experiment is obtained comparing n(k) in contrast to earlier findings based on comparing 
calculated and observed 3 He kinetic energies. 

PACS numbers: 67.60.-g, 61.12.Ex 



I. INTRODUCTION 

Liquid 3 Hc- 4 He mixtures are excellent examples of in- 
teracting Bose and Fermi liquids in nature. Their behav- 
ior is dominated by quantum-mechanical exchange effects 
which involve not only identical particles, as in the pure 
liquids, but also interactions between particles obeying 
different statistics. As a result, a rich array of macro- 
scopic properties are observed as the 3 He concentration 
is increased. Among these properties are a gradual sup- 
pression of the 4 He superfluidity and a finite solubility 
(~6.6%) of 3 He in 4 He at absolute zero temperature. 

Historically, investigationsi*2i2i4 of this binary mixture 
involved thermodynamical and hydrodynamical experi- 
ments. A central goal was to measure the phase diagram 
of the mixture and to investigate the elementary excita- 
tions as a function of 3 He concentration, x. In particular, 
the 3 He quasiparticles and the 4 He phonon-roton (p-r) 
spectral behavior have drawn much attention££iL& 

Fifty years ago, Landau and Pomeranchuk (L-P) pro- 
posed a model of the dilute helium mixtures at low tem- 
peratures as two interpenetrating fluids of 3 He and He 
atoms for which 3 He- 4 He interactions are retained and 
3 He- 3 He interactions are ignored. In this simple picture, 
the 3 He quasiparticle energy is given by es(k) = Hk 2 /2m* 
where m* is the 3 He effective mass in the presence of 4 He 
atoms. At low k, the 3 He excitation spectrum observed 
in neutron scattering^ agrees well with the L-P spectrum. 
However, for k > 1.5 A -1 , it departs from the theory and 
falls below the L-P curve. This observation and the ob- 



served small shift in the 4 He roton energji^, A, suggest a 
dynamical level-repulsion between es(k) and 64 (fc) in the 
vicinity of A. A smaller shift in the roton energy was also 
observed in Raman scattering 7 - for several mixtures con- 
centrations, up to x = 30%. As a possible explanation, 
Pitaevski: 9 proposed the existence of a small roton-like 
minimum in €3 (k). This suggestion found some theoreti- 
cal support^ but no experimental evidence. 

A fundamental parameter of interest in helium mix- 
tures is the energy-dependent quasiparticle effective 
mass, m*. Yorozu and collaborators^ reported high ac- 
curacy measurements of the 3 He effective mass in mix- 
tures and found an m* that strongly depends on pressure 
but varies little with x. Their result was later supported 
by calculations^ Recently, the effect of m* on the dy- 
namics of pure liquid 3 He has been investigated^ 

The first neutron scattering measurements at low mo- 
mentum transfer, Q, aimed at probing elementary exci- 
tations in helium mixtures were reported by Rowe et al& 
for a 5% mixture at saturated vapor pressure (SVP) and 
a temperature of 1.6 K and by Hilton et al& for several 
mixtures at SVP for temperatures in the range 0.6-1.5 K. 
In these pioniering investigations, the cryogenic capabil- 
ities were limited to temperatures close to, or above, the 
Fermi temperature, Tp . In this limit, thermal broaden- 
ing is important and can mask the quasiparticles energy 
band. Both Rowe et al. and Hilton et al. reported a 
small shift of the 4 He p-r curve but disagreed on the sign 
of the shift. 

Fak et al& were the first to investigate the excita- 
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tions at temperatures below Tp. Their measurement 
confirmed the small shift in the p-r spectrum observed 
by Rowe et al. and by Hilton et at. They found no 
evidence of a roton-like minimum in the 3 He spectrum. 
Their experimental results are consistent with the theo- 
retical predictions of Bhatl*^ and of Gotze et al.w^ 

Other complementary experimental techniques on he- 
lium mixtures include several studies of the elementary 
excitations by Raman scattering 7 ^ and a study of the 
static structure factor S(Q) by X-ray scattering^. We re- 
fer to the review articles by Ahlers 3 , Baym and Pethick 18 
and Glyde and Svensson-S and to the recent book by 
Dobbs 20 for further details on the subject. 

In recent years, inelastic neutron scattering at high 
Q has been widely used to measure the conden- 
sate fraction and atomic kinetic energy of pure liquid 
4 Hei 21 ' 22 i 23 i 24 i 25 ' 26 These high Q measurements are very 
challenging because of limitations imposed by Instrumen- 
tal Resolution (IR) and Final-State interactions (FS). A 
general review of the method and experiments on the 
VESUVIO neutron spectrometer is presented by An- 
dreani et al.SL In experiments involving 3 He, the large 
3 He neutron absorption cross-section reduces further the 
scattering intensity and thus the statistical precision 
of the experimental data. Nevertheless, a number of 
measurements at high Q have now been performed on 
pure liquid 2S i 29 i 30 i 31 ' 32 and solid 3 He^ 3 - In these measure- 
ments, sample cells having advantageous geometries with 
thin samples or samples at very low 3 He concentrations 
are generally used to minimize the effect of absorption. 

In contrast to the pure liquids, there are few reported 
neutron scattering measurements at high Q of the iso- 
topic helium mixtures. Most report values for the single- 
particle kinetic energies as a function of 3 He concentra- 
tion, x, in or close to, the normal phase 2 ^ 3 ^^ 1 .. The only 
neutron scattering measurement of no, to date, was re- 
ported by Wang and SokolM in a 10% mixture where 
they found an no = 18% ± 3. Several calculations of 
the Bose-Einstein condensate fraction (BEC), no, in he- 
lium mixtures 3 ^* 3 ^ 3 ^* 3 ^ 3 ^ predicted an enhancement of 
no- However, the enhancement predicted by these calcu- 
lations is significantly less than reported by Wang and 
Sokol. 

Theoretical treatments predict that the kinetic 
energy of both isotopes decreases with increasing 
concentrationi 35 i 39 i 40 i 41 Neutron scattering result a 29 ' 30 ! 31 
agree both quantitatively and qualitatively with theory 
for the behavior of if 4. In pure liquid 3 He, there is a 
reasonable agreement between recent experiments 29 ' 30,31 
and theory 4 ^ on the value of K3, although earlier ob- 
served values were lower than calculated values^ 3 . Maz- 
zanti et a/J-L have recently obtained good agreement be- 
tween theory and experiment by comparing a calculated 
dynamic structure factor J(Q,y) directly with the ob- 
served J(Q, y). 

In dilute helium mixtures, in contrast, there is a sig- 
nificant disagreement in K3. For example, the calculated 
value of K3 for a 6.6% mixture is about 19 K [Ref. [3|| but 



the measured value (^10-12 K) is found to be the same 
as in the pure liquid 3 He independent of x. Repeated re- 
finements of the calculation technique using methods as 
diverse as Monte Carlo diffusion calculations^ 5 - to path 
integral MC techniques^ 1 - have failed to resolve this dis- 
agreement. 

In this paper, we report on deep inelastic neutron 
scattering measurements aimed at determining no, if 4, 
and the 3 He momentum distribution in 3 He- 4 He mix- 
tures. The quantity observed is the dynamic structure 
factor (DSF) S(Q, ui). For a sufficiently large momentum 
transfer from the neutron to the struck atom, Q — > 00, 
the observed DSF reduces to the impulse approximation 
(IA) 2 ^, 

S IA (Q, u) = j dkn(k)S(cu -us.- k.w R ). (1) 

where oj r — HQ 2 /2m and v R = HQ/m are the free atom 
recoil frequency and velocity, respectively. In the IA, 
Sia(Q,w) depends only on a single 'y scaling' variable 
y = (u> — ll>r)/vr and is conveniently expressed as 



Jia(v) = v R Si A (Q,u) = J dkn(k)S(y - k Q ) 



(2) 



where fcg = k.j^r. JjA{y) is denoted the longitudinal 
momentum distribution. At finite Q, however, the struck 
atom does not recoil freely but rather interacts with its 
neighbors. These interactions introduce a Final-State 
(FS) b roadening function R(Q,y) in J(Q,y) (see, e.g. 

for a detailed account of FS effects). Including 
these final-state interactions (FS), the exact J(Q,y) at 
high Q is 



J(Q,y)= / dse- isv J(Q, s) = / dse iys JiA(s)R(Q,s) 

. (3) 

where J(Q, s) is the intermediate scattering function, 
Jia(s) is the Fourier transform (FT) of Jia{v) and is 
the one-body density matrix (OBDM) for displacements 
along Q, and R(Q, s) is the FT of the FS function 
R(Q,y). The reader is referred to Boronat et al~ and 
Mazzanti et alM> for a detailed discussion of S(Q,u>) for 
the isotopic helium mixtures. 

When there is a condensate, the fraction of atoms with 
zero momentum transfer, no, contributes an unbroadened 
peak to Jia(jj) and can be directly extracted from the 
data. A central goal is to find an appropriate model for 
the momentum distribution n(k) including a condensate. 
For quantum systems, n(k) is also not a Gaussian. In 
earlier work 2 ^ -, we found that n(k) in liquid He and 
in solid 4 He differs markedly from a Gaussian. 

This report is laid out as follows. The experimental 
method is presented in section[Hj followed by a discussion 
on the data reduction and analysis methods in section 
IIIII The results are presented in section llVl followed by 
a discussion and conclusion in section Ivl 
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FIG. 1: MARI Data showing J(Q, y) for pure liquid 4 He in the nor- 
mal (open circles) and superfluid phase (closed circles). A higher 
peak arising from the Bose-Einstein condensate fraction is observed 
in the superfluid phase. The instrument resolution function is over- 
plotted (dotted line) for comparison. 



A special sample cell having a slab geometry, of the 
type described by Sokol et al£& and by Fak et al&, was 
designed. The sample cell was then placed in the beam 
and cooled using a 3 He sorption cryostat. Two separate 
experiments were performed on MARI under identical ex- 
perimental conditions. For each set of measurements, a 
separate background measurement of the empty cell was 
taken, and the two were found to be very similar, as ex- 
pected. Data was collected at temperatures of T = 0.4 K, 
T = 1.3 K and T = 2.5 K which were measured using Ge 
temperature sensors located at the bottom and top of the 
sample cell and connected to a Neocera temperature con- 
troller. A gas handling system at room temperature was 
used to prepare mixtures of 3 He gas to yield liquid 3 He 
concentrations of x — 0, 5, 10, 15, and 20 %. The pres- 
sure above the liquid was maintained at saturated vapour 
pressure, and was monitored with a pressure transducer 
by keeping a capillary line open between the gas handling 
system and the cell. 



B. Data Reduction 



II. EXPERIMENTAL DETAILS 

A. Experiment 

The experiment was performed using the MARI time- 
of-flight (TOF) spectrometer at the ISIS spallation neu- 
tron source at the Rutherford Appleton Laboratory, 
United Kingdom. MARI is a direct geometry chopper 
spectrometer in which the time of arrival of a neutron 
in the detector, measured from when the neutrons leave 
the moderator, determines its energy loss or gain after 
scattering from the sample. The momentum transfer de- 
pends on both the TOF of the neutron and its scattering 
angle. More than 900 3 He gas detectors provide a cover- 
age of scattering angles between 3° and 135° in steps of 
0.43°. A large range of momentum and energy transfer 
can therefore be observed simultaneously. 

An incident neutron energy of 765 meV was selected 
allowing wavevector transfers up to Q = 30 A" 1 and en- 
ergy transfers up to E = 700 meV to be investigated. 
Measurements of J(Q,y) were made for wavevectors 20 
< Q < 29 A -1 . The low energy resolution Fermi chop- 
per was selected to increase count rate in the presence 
of the highly absorbing 3 He. This has an energy resolu- 
tion of approximately 25 meV. The resolution function 
(in y) was determined at each Q as follows. The intrin- 
sic J(Q,y) of pure liquid 4 He is accurately known from 
previous measurements on MARI using a high resolution 
setting. 26 Using this J{Q, y) as input, the present broader 
instrument resolution can be determined directly by re- 
producing the observed resolution broadened J(Q,y) of 
pure He {x = 0). The resolution determined in this way 
at Q = 27 A -1 is shown in Fig. [T] 



Standard procedures were employed to convert the 
raw neutron scattering data from TOF and scattered 
intensity to energy transfer and the dynamic structure 
factor, S(Q,w). A summary of these conversions has 
been given by Andersen et al.^L The data was then con- 
verted to the y scaling energy transfer variable and to 
J{QiV) = vr.S(Q,u>). The data were analyzed in detail 
for wave vectors 26 < Q < 29 A -1 . 

The background arising from the empty cell was mea- 
sured separately and found to be a smooth function of 
scattering angle and TOF with no microstructure. At a 
given Q value, the measured background was well repre- 
sented by a smooth quadratic function in y, ay 2 + by + c, 
where a, b and c were determined by a least square fit 
to the measured background. An example of this back- 
ground is shown in Fig. |2]as a dotted line. This measured 
background function was subtracted from the scattering 
from the helium plus cell to obtain the net measured 
scattering intensity from liquid helium. 

The analysis consists of representing the net J(Q, y) as 
a sum of models for the scattering 4 He and 3 He. A linear 
background of variable height and slope was also included 
to allow for any error in the background determination 
or shielding of the background by the sample. This lin- 
ear background was always very small. Additional error 
in model parameters for 4 He and 3 He arising from incor- 
porating the linear background is included in the quoted 
parameter errors below. The data was y scaled so that 
the 4 He peak is centered at y = using the 4 He mass. 
This means that the 3 He peak component of the model 
function must be multiplied by M — m^jm^ to regain 
the correct magnitude. An example of the total scat- 
tered intensity y scaled and of models 4 He and 3 He plus 
the measured background is shown in Fig. [21 
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FIG. 2: Total observed scattering intensity (solid circles) showing 
J(Q> y) f° r a 3 He- 4 He mixture at 10% 3 He concentration and T = 
2.5 K plus background. Dotted line is a separate measurement of 
the background. The data is presented with the 4 He peak centered 
at j/4 = 0. The error bars are the size of the circles. 



III. DATA ANALYSIS 

Our goal is to determine the condensate fraction and 
to get as much information as possible about the He 
and 3 He momentum distributions in the mixture. The 
data is good enough to determine at most two free fit- 
ting parameters for each peak. To describe the 4 He peak, 
we follow the method used for the pure liquid 4 He. In 
superfluid 4 He, we use the convolution approach (CA) 
in which a FS broadening function R(Q, y) is convoluted 
with a model Jja(h) as in (JSJ. In normal 4 He, we use the 
additive approach (AA), described below, and the CA. 
The condensate fraction and the width (~ kinetic energy) 
of the peak are the free parameters. The FS function and 
the shape of n(k) are assumed to be the same as in pure 
4 He. The comparison between normal and superfluid liq- 
uid 4 He in Fig. 2] shows that the peak height is higher 
at y = for the superfluid phase due the presence of 
a condensate and a left-right asymmetry introduced by 
the condensate term noR(Q,y) in J(Q,y). To describe 
the 3 He peak, we construct a model of a Fermi momen- 
tum distribution n(k) and compare the corresponding 
resolution broadened Jia (y) in (0) with the data (no FS 
effects). The parameters in the 3 He n(k) are the step 
height at the Fermi surface and the length of the high 
momentum tail. 

There is a small overlap of the net 3 He and 4 He peaks. 
Fig. [21 shows that the degree of overlap is very small 
indeed. The area of overlap of the models is typically 
1-2% of the total area under the two models. We found 
that the parameters in the model fit to one peak were 
independent of the model used for the second peak within 
the error quoted. This was tested by using various models 
to represent the second peak. 



A. 4 He response 

1. Convolution Approach 

The convolution approach (CA) is most useful when 
the momentum distribution n(k), and thus JiA{y) in J2J, 
are narrow in y relative to R{Q, y). Equally, it is useful if 
n(k) contains a component, such as a condensate, that is 
narrow relative R(Q, y). In this case, the observed width 
of the narrow component in J(Q,y) is set by R(Q,y). 
In pure liquid 4 He, the FWHM of R(Q, y) is approxi- 
mately I A -1 . This is broad compared to a condensate 
component but narrower than the remainder of n(k). 

To determine the condensate fraction, it is convenient 
to separate the state that is macroscopically occupied 
from the regular uncondensed states. In a uniform liquid, 
the natural orbitals are plane wave, momentum states. 
The orbital containing the condensate is the k = state. 
This state contributes a term nn<5(k) to n(k) where no = 
Nq/N is the fraction of particles in the condensate. The 
regular, uncondensed states are the k ^ states and 
we denote their contribution to n(k) as n*(k). In an 
interacting Bose liquid, Bosons can scatter into and out 
of the condensate from the k ^ states. This leads to 
a coupling between the condensate and the k ^ states 
and a terr^Z^ 
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in n(k). The coupling is strongest for the low k states 
and the expression in the square bracket of (@J is derived 
for and valid for the low k limit. We have multiplied it 
by a Gaussian to cut off no/ (k) at higher k as done for 
pure2& 4 He with k c fixed at k c = 0.5 A -1 . This term 
is highly localized around k = and when Final State 
broadening is included no/(k) cannot be distinguished 
experimentally from no<5(k). The model n(k) we use for 
the 4 He component therefore has three terms2i, 



n(k) = no[<J(k) + /(k)]+Am*(k) 



(5) 



Ai is a constant chosen by normalization J cZkn(k) = 1. 

To implement this model, we Fourier transform n(k) 
to obtain the corresponding One-Body-Density-Matrix 
(OBDM), n(r)=f dke lk r n(k). J IA (s) is the OBDM for 
displacements r = sQ parallel to the scattering wavevec- 
tor Q. Jia(s) corresponding to (J3J) is obtained by Fourier 
transforming JSJ to obtain n(r) and using Jia(s) = n(s) 
giving, 



Jia(s) = n(s) = n [l + f(s)] + Am*{s) 



(6) 



The term [I + f(s)] is long range in s. The n*(s) is 
short range in s and we represent it by2& 



n*(s) = exp 
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where 0,%, and ae are parameters (cumulants). This 
representation is useful when n*(s) is at least approxi- 
mately a Gaussian, plus small corrections. In going from 
iJSJ to |JBJ, we actually need to transform /(k) only, and 
once only since /(k) has no free parameters. Given /(s), 
our model may be viewed as a model for Jja{s) in 

The model Jia(s) is multiplied by the FS function 
R{Q, s) and the product is Fourier transformed to obtain 
J(Q, y) as in 10. In the present fit to the 4 He peak, we 
used the pure 4 He FS function at all 3 He concentrations, 



R{Q, s) = exp 



0: 
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(8) 



with parameters (3% = 0,3/ XQ, (3$ = a^/iXQ) 3 and 
= &64/(AQ) 2 with (13,052 and set at their pure 
4 He values^ We also set 0:4 and 6tQ in n*(s) at their 
pure 4 He values^ The J(Q,y) therefore has only two 
free parameters, ai and uq. 

The above procedures for determining the shape of 
the momentum distribution of pure liquid 4 He, and ex- 
tracting the condensate fraction are well established, 
and have been the subject of intense activity over the 
years, 21 . 26 ! 47 . 52 



2. Additive Approach 

When the momentum distribution is broad, FS effects 
are relatively less important. In this case, the exponential 
in R(Q, s) of JSJl can be expanded and FS effects retained 
as additive corrections to the IA. Similarly, if the devia- 
tions of n*(s) from a Gaussian are not too large, n*(s) 
in can be expanded and deviations from a Gaussian 
n*(s) retained as additive corrections. This expansion 
of n*(s) is the same as the Gauss- Hermite expansion in- 
troduced by Sears* 2 ^ These expansions lead to the Addi- 
tive Approach (AA), in which J(Q,s) is represented as 
a Gaussian IA, Jja{s), plus corrections for deviations of 
n*(s) from a Gaussian and for FS effects a a 2 ^ 53 , 



J(Q, s) = J IA (s) + Ji(Q, s) + J 2 {Q, s) 



(9) 



where 



Jia(s) = (27ra 2 ) 
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M s ) = y^3S 3 JlA(s) 
J2(s) = ^4S 4 J/a(s) 



(10) 
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FIG. 3: Schematic diagram of the Fermi model momentum dis- 
tribution used in this analysis. The parameters are defined in the 
text. 



B. 3 He Response 

The 3 He momentum distribution is not well approx- 
imated by a Gaussian. Thus we do not expect the fit 
of a Gaussian or the A A to the 3 He peak to reveal the 
physics well. To proceed, we constructed a simple model 
of the 3 He n(k) consisting of a step of height Z at the 
Fermi surface plus an exponential high momentum tail. 
We calculate the Jia{v) for this n(k) using J2J and fit 
the convoluted Jia(u) directly to the data. 

Specifically, the model n(k) is, 



n(k) 



H 

V F 



-r- < 1 
r- > 1 



(11) 



where kp is the Fermi momentum and Vp — ^§-kp the 
volume of the Fermi sphere. Assuming that the volume 
occupied by 3 He in the mixtures exceeds that occupied 
by 4 He by about 28%, as suggested by the dielectric 
measurements of Edwards et al^, the Uf can be ap- 
proximated using the expression kp = (37r 2 px) 1 ^ where 
p is the average density of the mixture at SVP and x 
the 3 He concentration. The kp values are listed in Ta- 
ble [I] The model has three parameters, H, Hp and 
77 where 77 determines the length of the high momentum 
tail and Z = H - Hp (See Fig. 0). Normalizat ion of 



and P3 = CI3/XQ and pi — 0:4 + a^/(XQ) 2 . The 6:4 is the 
leading deviation from a Gaussian and 03 and 04 are FS 
terms. In fitting this function to the 4 He peak in normal 
mixtures, we fixed 03 and 0:4 at their pure 4 He values 
(S4 = 0). The J(Q, y) then has only one free parameter, 
CI- 



TABLE I: Fermi momentum kp(x) = (37r 2 px) 1 / 3 as a function of 
the 3 He concentration x where p is the average mixture density at 
SVP. 
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n(k), J dkn(k) = 1, gives 

3H F ( 2 2 \ 

H+ + _ + l (12 

rj \r) 2 7] J 

which can be used to eliminate one parameter, say H. 
Similar models of n(k) are discussed by Carlson et a/4* 
and used by Azuah et al£& and Mazzanti et aim* for pure 
liquid 3 He. The IA given by corresponding to (|llfl is, 

7 , , fw(*(w 2 ) + ^(i+±)) y'< 1 

(13) 

where y' = y/kp. 

If the liquid temperature T is near or above the Fermi 
temperature, Tf = 2Tr f» ks kp , there may be some thermal 
broadening of the Fermi liquid n(k). 20 This is certainly 
the case for a Fermi gas for which n(k) is a step func- 
tion (Z = 1) at T = K. However, in the liquid model 
above, interaction has already reduced the magnitude of 
the step at kp (Z < 1). Thus thermal broadening may be 
relatively less important in a strongly interacting liquid. 
This is discussed more fully in section llVl 

To estimate Tf, we note that the effective mass m* 
depends strongly on the pressure but little on the 3 He 
concentration^ We take m* — 2.3m as found by theoret- 
ical calculationsi2^£, independent of concentration. In 
the range x = 5-20 %, the corresponding Tp is 0.3-0.8 K. 
This is significantly less than the normal liquid temper- 
ature T = 2.5 K and comparable to the superfluid tem- 
perature T — 0.4 K used here. Our n(k) extracted from 
data will represent n(k) at these temperatures. However, 
we did not observe any difference between n(k) at T = 
0.4 K and 2.5 K. This suggests that the rounding of the 
small step Z at kp by FS effects and the instrument res- 
olution function is comparable or greater than thermal 
broadening. 



IV. RESULTS 

In this section, we present our neutron scattering data 
and the results for the 4 He atom kinetic energy, the 3 He 
momentum distribution and the 4 He condensate fraction 
obtained from analyzing the data. Since the absorption 
cross section of 3 He is so large, the net scattered neutron 
intensity from 3 He- 4 He mixtures is weak. For this reason, 
a broad energy resolution setting was used on the MARI 
spectrometer, as shown in Fig. ^ Also, the statistical 
precision of the data is not high. As a result, we were 
able to determine uniquely only one or two parameters 
in model fits to the data. Specifically, we found that 
the parameter no that provides the condensate fraction 
and the parameter cxi that sets the kinetic energy, K4 = 
3h 2 ct2/2m, were correlated. For this reason, we begin 
with normal 4 He where uq = 0. 



A. Normal phase 

1. 4 He kinetic energy 

The kinetic energy K4 of 4 He in the normal liquid at 
T = 2.5 K is shown in Fig. 0|as a function of 3 He concen- 
tration. K4 is determined from the 4 He peak using two 
fitting procedures, the convolution approach (CA) and 
the additive approach (AA). In the CA, we use the Final 
State function that was determined^ previously for pure 
4 He unchanged at all 3 He concentrations. The only free 
parameter in the present fit was 612 ■ Similarly, in the AA 
the FS parameters were set at their pure 4 He values 2 ^ 
and only 0,2 was free. The parameters held fixed in the 
CA are summarized in Table ITT1 

Fig. 21 shows that the CA and AA procedures give 
identical kinetic energies for pure 4 He (x = 0). This 
value agrees with our previous^ 4 He determination 
K4 =16.3±0.3 K and with the value found by Sen- 
esi et al.m^ As the 3 He concentration is increased, the 
two methods continue to agree within error but the CA 
method gives a marginally lower K4. The K4 decreases 
only very little with increasing x. Azuah et alm& and 
Senesi et alM find that K4 decreases more sharply with 
increasing x. A decrease is anticipated since the density 
decreases as x increases and in pure fluids K4 decreases 
as density decreases. When we fit a Gaussian to the 4 He 
peak, we find a K4 that is independent of x. We return to 
this point in section TlV CI and compare with calculations 
in section IVI 

The K4 values in Fig. 0]are averages obtained from fits 
to data at specific Q values. The variation of K4 with Q 
is shown in Fig. [S] The aim is to display the statistical 
precision of the data. 



2. 3 He momentum distribution 

In this section, our goal is to learn as much as possible 
about the 3 He momentum distribution. We do this by fit- 
ting the JiA(y) (|13|) obtained by substituting the model 
n(k) given by l|ll|) and shown in Fig. |3J into 10) • The 
model n(k) has two parameters, Hp giving the height 
of the high momentum tail at kp and 77 which sets the 
length of the tail. H is determined by normalization. 



TABLE II: Fitting variables used in the CA. In the normal phase, 
only Q2 is determined from a fit to data and in the superfluid phase 
only no is determined. The values of fc c , j3 n and a„ for n > 3 were 
all kept fixed at their values in the pure 4 He case, as evaluated by 
the precise measurements of Clyde et al.ii 
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FIG. 4: 4 He kinetic energy [K 4 ] at T = 2.5 K as a function of 3 He 
concentration, x. The closed squares are the 4 He kinetic energies 
obtained from the Additive Approach (AA) and the open circles 
obtained from the Convolution Approach (CA)™ 



The data in 3 He peak region was precise enough to de- 
termine one parameter well with an estimate of a second 
parameter. Since the tail of n(k) is of specific interest, we 
also set the tail parameter at specific values e.g. no tail 
(H = 1), tail that correponds to the calculated n(k) 39 
(r) = 0.8) and fitted for Hp to see how good a fit could 
be obtained. We subsequently obtained best fit values 
of 77. Once, 77 is determined, the model kinetic energy is 
obtained as K3 — J dkfc 2 n(k) giving, 



3ft 2 fc| 
2m3 5 



H 



V V 



4! 



12 

1^2 



4 
V 



r 



(14) 

where ^— = 8.08 K A 2 . For typical values of 77, more 
than 95% of K3 arises from the tail in n(k). Also, in 
(|14fl . K3 is very sensitive to 77. Thus, unless 77 can be 
accurately determined, K3 is not a good single parameter 
to characterize the data. 

Fig. shows fits to the 3 He peak with no high mo- 
mentum tail (Z = H = I, Hp = 0) and with the tail 
parameters set at 77 = 0.8, the value obtained in a fit to 
the high k portion of the tail calculated by Boronat et 
al. Clearly a Fermi step function with no tail cannot 
reproduce the data. In contrast, an n(k) with a high en- 
ergy tail that reproduces the calculated tail (see Fig. EjJ 
fits the data. In this sense, the calculated n(k) and the 
data are entirely consistent. There is no disagreement 
between theory and experiment. 

We determined an observed value of 77. The best fit 
value of 77 obtained varied somewhat with the data set 
(e.g Q value) considered reflecting the statistical pre- 
cision of the data. Fig. shows fits for 77 = 0.76 
and 0.94 which provide equally good fits at different 
Q's. The data is precise enough to determine that 77 
lies in the range 0.7 < 77 < 1.0, i.e 77 = 0.85±0.15 and 
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FIG. 5: Top: Additive Approach (AA) fit to data in the normal 
phase (T = 2.5 K) at several 3 He concentrations showing the Q 
dependence of K4. Bottom: Convolution Approach (CA) fit to 
data in the normal phase. 



Z = H — Hp = 0.05 ± o of . This correspond to a very 
strongly interacting Fermi liquid. Fig. [S] shows a com- 
parison between our observed n(k) for a 5% mixture at 
T = 2.5 K and a calculated^ n(k) at 6.6%. 



B. Superfluid phase 

As noted, we found that the parameter cb and the 
condensate fraction no were correlated. For example, in- 
creasing no and decreasing both lead to a narrower 
peak that cannot be distinguished within the present pre- 
cision of the data. In pure 4 He, a.^ was found to be the 
same in normal (T = 2.5 K) and superfluid (T = 0.5 
K) 4 He. Thus in the mixtures we set ct2 in the super- 
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FIG. 6: Fits to the 3 He peak (solid lines) in the normal helium mixture data at x = 5% and Q = 29 A -1 . Data are open 
circles with error bars. The error bars are the size of the circles except where shown. The left side shows a fit assuming 
the 3 He rt(k) is a Fermi step function (i.e n(k) has no high momentum tail (H =1 )) a nd the right side shows a fit to data 
with the model n(k) shown in Fig. |J]that reproduces the calculated tail of Ref. I.'>!l . The fit with no tail (left) is poor 
while the fit incorporating the calculated tail (right) is good. 




FIG. 7: Fits to the 3 He peak (solid lines) in the normal helium mixture data (open circles) at x =5% and Q = 28.5 
A" 1 . The error bars are the size of the circles except where shown. The graphs show that it is possible to get good fits 
to the same data with two different sets o f pa rameters. This illustrates the precision in which the n parameter can be 
determined. The calculated tail from Ref. l3fll corresponds to r\ = 0.8. 



fluid phase at the normal phase value determined above 
at each x. 

1. Condensate fraction 

To determine no (x) , we first perform fits of the C A to 
data for each Q at a given concentration x. The average 
value over all Q gives the corresponding no (x) . The varia- 
tion of the condensate fraction with Q is shown in Fig. [5] 
The variation reflects the statistical precision of the data. 
Fig. 1101 shows the resulting x dependence. The helium 
mixtures data at x = 20% was taken at a slightly higher 
temperature (T = 1.3 K) than the lower concentrations 



data. To be able to make a consistent comparison of the 
results on n Q , we made a temperature correction^ to n 
at x = 20% using n (T) = n [l- (^) 7 ] where T\ - 2.17 
and 7 = 5.5 to get its equivalent value at T = 0.4 K. The 
corrected n Q is shown in Fig. ^|along with the results at 
lower 3 He concentrations. We find that no (a;) increases 
slightly above the pure 4 He value with increasing x. Fig. 
1101 represents our final results for no which are discussed 
in section fvl 
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FIG. 8: Comparison of the observed 3 He atom momentum distri- 
bution n(k) (solid line) at x =5% with the calculated n(k) from 
Ref. 15s! at x = 6.6% (dotted line). 



2. 



i He kinetic energy 



The 4 He kinetic energy in the superfluid phase can be 
obtained using the form of the momentum distribution 
in ijjjj and assuming that the parameter a.i is the same in 
both normal and superfluid phases. Normalizing the He 
model momentum distribution, as discussed above, leads 
to n [l + If] + At = 1 with I f ~ 0.25. This result can 
be obtained from © as n(s — 0) where If = f(s = 0) 
and n*(s = 0) = 1. Since n Q is a function of x, the 
normalizing constant Ai is also a function of x, Ai(x) = 
1 — 1.25tt,o(.t). The kinetic energy K4 is proportional 
to the second moment of n(k). Only n*(k) contributes 
significantly to the second moment so that from JSJ and 
with a2 the same in the normal and superfluid phases we 
have K^(S) = A%(x)K4(N). In this model, the kinetic 
energy in the superfluid phase (S) is reduced below that 
in the normal phase (N) entirely by BEC. We found the 
kinetic energy to be about 10% lower in the superfluid 
phase than in the normal phase increasing somewhat with 
x as shown in Fig. 1111 



3. 3 He momentum distribution 

Fig. 1121 shows two fits to the 3 He peak in the super- 
fluid mixture at T = 0.4 K. In one fit, the tail parameter 
of the model n(k) is set at 77 = 0.8, the value that re- 
produces the calculated tail. The same model (77 = 0.8) 
was compared with data for T = 2.5 K in Fig. [fj] The 
agreement with the data for 77 = 0.8 at T = 0.4 K and 
2.5 K is indistinguishable. This shows firstly that there is 
no observable temperature dependence of the 3 He peak, 
between T = 0.4 K and T — 2.5 K which corresponds to 
T ~ Tp and T ~ 6Tp, respectively. 
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FIG. 9: The 4 He condensate fraction (closed squares), no, as a 
function of momentum transfer Q. The error bars represent statis- 
tical errors obtained from fits to data. The solid lines are guides 
to the eye. 



Somewhat surprisingly, there is no observable temper- 
ature dependence of the 3 He momentum distribution al- 
though the two temperatures are significantly different. 
We return to this point in the discussion section. Sec- 
ondly, there is again no disagreement between theory (for 
T = 0K) and experiment. 

As at T = 2.5 K, a range of r\ values provided best fits 
to data depending upon the Q value considered. Fig. IT21 
shows an example at Q = 27.5 A -1 where the best fit is 
obtained for r\ = 1.0. At T — 0.4 K, we find r\ lies in the 
range 0.7< r] < 1.0 or 77 = 0.85±0.15 as found at T = 2.5 
K. 



C. Alternate methods of analysis 

In this section, we test the sensitivity of the results to 
the method of analysis by analyzing the data in different 
ways. 
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FIG. 10: The 4 He condensate fraction (closed squares), no, as a 
function of 3 He concentration. The error bars are standard devia- 
tions obtained from a linear least-squares fit to the Q dependence 
of no- The condensate fraction increases with increasing x from 
7.25± 0.75% at x = to 11±3 at x = 15-20%. The dashed line is 
a guide to the eye. 
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FIG. 11: Kinetic energy of 4 He atoms, K4, as a function of 
concentration. The superfluid K4 were obtained from the 
normal K4 assuming the condensate does not contribute to 
the kinetic energy. 



1. Kinetic energy assuming a Gaussian 

A straightforward method of analysis is to fit Gaussian 
functions to both the 3 He and 4 He peaks as employed by 
Azuah et al.w& This method, when applied to our data 
yields a 3 He kinetic energies that agree with those previ- 
ously obtaine d 29 i 31 where a Gaussian or a Gaussian-like 
fit was made. Senesi et alM^L also included the leading 
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FIG. 12: Fits to the superfluid helium mixture data at a 3 He 
concentration x = 5%. The 4 He peak is fitted using the CA method 
and the 3 He peak assuming a Fermi n(k) having a variable tail 
length n. The open circles are data points and the solid lines are 
fits to data. The error bars are the size of the circles except where 
shown. 



term beyond a Gaussian in a Gauss-Hermite expansion 
of J(y) to analyze their data. However, they found that 
the inclusion of higher order terms did not result in sig- 
nificant improvements in the fits. Both Azuah et al. and 
Senesi et al. found K4 values that decrease with increas- 
ing 3 He concentration. 

In Fig. El we show the K4 values obtained by fitting a 
Gaussian J(Q, y) to our data compared with the reported 
values by Azuah et alm& and by Senesi et al.m^ Firstly, 
for pure 4 He (x = 0) we find a K4 that lies approximately 
10% below the CA values. A smaller K4 is expected for 
a Gaussian. In the CA, FS effects cut off the wings of 
the peak and allow a somewhat broader ?i(k) to fit the 
data. Senesi et al. took data at Q ~ 100 A -1 where 
FS effects are probably negligible. Thus we expect their 
values to agree with our CA values. Most importantly, 
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FIG. 13: Observed K4 (open circles) assuming a Gaussian fit to 
the 4 He J(Q,y). The measured values by Azuah et al. (open 
triangles) 29 and Senesi et al. (open squares)^™ are shown for com- 
parison. The dashed lines are guides to the eye. 



FIG. 14: Condensate fraction, no, obtained from present data 
assuming the kinetic energies measured by Azuah et al. and Senesi 
et al£L. The no decreases with increasing x. This is illustrates the 
correlation between no and K4. 



using a Gaussian fit we find a K4 that is independent of 
x. Thus while our CA, the Azuah et al. and the Senesi et 
al. K4 agree for pure 4 He value (x = 0), the data itself 
appears to disagree on the x dependence of K 4 . Most 
theoretical calculations also find that K4 decreases with 
increasing x as shown in Fig. 115b . 

2. Condensate assuming concentration dependent kinetic 
energies 

Since there is some correlation between the c\i param- 
eter (if 4) and no, we determined no using the values of 
CY2 corresponding to the measured 4 He kinetic energies 
by Senesi et al.~ The resulting values of no are shown 
in Fig. 1141 We find no decreases with increasing x, when 
the x dependent K4 values of Senesi et al. and Azuah et 
al. are used. This result is expected since a narrower 4 He 
peak is obtained by either increasing no or decreasing 0.2- 
If c<2 (~ K4) decreases, a smaller no is needed to obtain 
a good fit. Thus a K4 that decreases significantly with x 
implies an no that also decreases with x in disagreement 
with theory. 

V. DISCUSSION AND CONCLUSION 

A. 4 He response 

In pure normal 4 He, we find a kinetic energy [K4 = 
16.1 ± 0.3 K) which agrees well with previous neu- 
tron scattering data™ (K4 = 16.3 ± 0.3 K ) and with 
calculations^ (K4 = 15.41 K). In normal helium mix- 
tures, we find a K4 that decreases somewhat with in- 



creasing 3 He concentration. This suggests that the lo- 
cal environment of individual 4 He atoms remains largely 
unchanged^ with the addition of 3 He atoms at low x. 
This finding represents the chief difference between our 
results and early theoretical treatments" and previ- 
ous experiments 29,31 which find that K4 decreases signif- 
icantly with increasing x. Theoretical and experimental 
values of K4 are compared in Fig. 115b . The present K4 
and a more recent Path Integral Monte Carlo (PIMC) 
calculation^ of K4 which show a substantial agreement 
on the x dependence are compared separately in Fig. 115b . 
Clearly, there remain some theoretical and experimental 
differences to be resolved. The present K4 in the super- 
fluid phase was obtained from the normal phase value 
assuming that the condensate does not contribute to the 
kinetic energy. 

To determine the condensate fraction, no, in helium 
mixtures, we introduced a model n(k) which has a con- 
densate term and fitted the model to the neutron scat- 
tering data. In the model, the second moment of n(k) for 
the finite k states was assumed to be the same as in the 
normal phase, as we found previously 2 ^ for pure liquid 
4 He. The higher moments and the Final State function 
were assumed to be the same as in pure superfluid he- 
lium. 

In mixtures we find a small increase in no with in- 
creasing 3 He concentration, to no = 11 ± 3% at x = 
15-20 %. Fig. 1161 shows our observed no as a function 
of x compared theoretical values. Our results agree with 
the theoretical values within experimental error. Wang 
and Sokol 34 measured no at one concentration reporting 
a value of no = 18% at x = 10%. This is higher than 
the present measured and calculated values. The con- 
centration dependence of no may be roughly estimated 
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FIG. 15: (a): 4 He kinetic energy (K4) as a function of 3 He concentration, theory and experiment. Calculations are VMG™ (solid 
circles), HNC using a Jastrow function™ (triangles), HNC including triplet correlations^? (solid squares). Experiments are Azuah et alm& 
(up-pointing triangles), Wang and Soko 134 (down-pointing triangle), Senesi et aldil (open squares) and present at T = 2.5 K (open circles). 
Dashed lines are fits to experiments, (b): Present observed K4 at T = 2.5 K (open circles) as a function of 3 He concentration compared 
with a PIMC calculation™ of K4 at T = 2.0 K (solid stars). Dashed lines are guides to eye. 



assuming that the effect of the 3 He is simply to change 
the volume available to the 4 He (average density approx- 
imation) and using the density dependence of no in pure 
4 Hei— 1 Essentially, the molar volume of the mixture at 
low concentrations is approximately— v ~ V4 (1 + ax) 
where V4 is the 4 He molar volume at SVP and a the ex- 
cess of volume occupied by 3 He (~ 0.28). This gives a 
mixture density which decreases with x. On the other 
hand, the condensate fraction, no, in pure liquid He 
increases22i££ with decreasing density. This yields an no 
which increases very slightly with x, from no ~ 7% (x — 
0) to n ~ 9% (x = 20%). This simple calculation sug- 
gests that no should increase little with concentration, in 
qualitative agreement with our experiment. 
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We find that a model 3 He momentum distribution n(k) 
with a small step, Z, at the Fermi surface (k — kp) and 
a substantial tail at higher k characteristic of a strongly 
interacting Fermi liquid reproduces the observed 3 He re- 
coil peak well. For a gas of non-interacting Fermions at 
T = K, n(k) is a simple step function (Z —1); n(k) = 1 
for k < kp and n(k) = for k > kp. This case is repro- 
duced by our model n(k) in Eq. (|ll(l for parameters H = 
1, Hp = (Z = H - Hp = 1, no tail in n(k)). When 
Fermions interact, some Fermions move from states be- 
low kp to states above kp to minimize the total energy, 
potential plus kinetic energy. The step at kp is reduced 
Z < 1 and n(k) develops a tail at higher k. This is repro- 
duced in Eq. for parameters H — Hp = Z < 1 and 
77 > 0. The stronger the interaction, the smaller is Z and 



FIG. 16: The condensate fraction (no) as a function of 3 He con- 
centration, theory and experiment. Calculations are the open sym- 
bols; VMG™ (circles), HNC (Lennard- Jones )Sm (down-pointing tri- 
angles), HNC (triplet correlations )22 (squares), HNC (Jastrow)ffi 
(up-pointing triangles), HNC (Aziz potential) 3 ?, (right-pointing tri- 
angles). Measurements are Wang and Soko! 34 (solid circle) and 
present (solid squares). 



the longer is the tail (smaller 77). In pure liquid 3 He a 
recent diffusion Monte Carlo (DMC) calculation 4 ^ finds 
Z = 0.236 and a tail in n(k) that is well represented by 
an exponential (times algebraic terms) in k. This n(k) 
agrees well with earlier DMC results 

Pure liquid 3 He is regarded as a strongly interacting 
Fermi liquid, more strongly interacting than nuclear mat- 
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ter, and Z ~ 0.25 represents a large reduction of the step 
height at kp arising from interaction. As noted above, 
Boronat et al£^ calculate an n(k) obtaining Z = 0.09 and 
a tail that is also well fitted with an exponential with 77 = 
0.8 for a 3 He- 4 He mixture of x — 6.6 % 3 He. Apparently, 
3 He in a dilute mixture is more strongly interacting than 
in pure 3 He. Our observed, finite temperature n(k) has 
Z = 0.05 ± Ej^ 4 and r] = 0.85 ± 0.15 at x = 5 %. It agrees 
with the calculated n{k) for T = K. Z = 0.05 ± g;Jf 
represents a very small step at kp . 

We have not included Final State (FS) effects in our 
analysis of the 3 He recoil peak. The impact of FS ef- 
fects depends critically on the width and shape of the 
intrinsic Jia{d) arising from n(k). If Jia{v) is narrow 
relative to the FS broadening function R(Q, y) or has a 
narrow component or sharp feature, then the FS effects 
serve to broaden the narrow component to the width of 
the FS function. The classic example is the condensate 
component in superfluid 4 He. However, if Jja(]j) is broad 
relative to the FS function, then the FS function has little 
further broadening impact. For example, the second mo- 
ment of the FS function is zero. In this event the FS ef- 
fects can be well represented by a moment expansion 4 ^^ 
in which the leading term is the third moment. This term 
serves to introduce an asymmetry into J(Q,y) but will 
not modify its width or whether J(Q, y) has significant 
tails at larger y arising from tails in n(k). 

The FWHM of R(Q, y) at Q w 25 A" 1 is in the range 
0.7-1.0 A -1 in pure liquid 3 He, in 3 He- 4 He mixtures and 
in pure liquid 4 He. 26 i 44 i 49 The FWHM of Ji A (y) of 3 He 
in 3 He- 4 He mixtures is approximately 2.0 A -1 , largest 
in dilute 3 He mixtures. We expect FS effects firstly to 
spread out the discontinuity in Jia(jj) arising from the 
step of height Z in n(k). It will be spread out over a 
width y ~ 0.7 A -1 (i.e. of order of 2kp). The present 
instrument resolution function, shown in Fig. ^ which is 
broader than the FS function at Q ~ 25 A -1 , will further 
spread out the discontinuity. Parenthetically, it is prob- 
ably largely because of the resolution and FS broadening 
that we do not observe any apparent broadening of n(k) 
arising from temperature. Otherwise, the FS function is 
narrow relative to the width of Jia(d) arising from n(k). 
Thus we do not expect FS effects to change the overall 
width of J(Q,y) (i.e. the apparent width of n(k) and 
whether it has tails or not). However, including FS ef- 
fects should improve the shape of J(Q, y) as it has been 
done^ in pure liquid 3 He. Thus we do not expect FS 
effect to alter our basic conclusions about n(k). 

We note that one cannot use the FS function for pure 
3 He for mixtures. Essentially, in dilute 3 He mixtures the 
Fermi statistics for only the 3 He atoms should be in- 
cluded rather than for all atoms as in pure 3 He. We plan 
to incorporate FS effects in a future publication. 

We have strictly used a model n(k) valid for T = K 
since the step Z at kp is not rounded by thermal effects 
to fit data taken at T = 0.4 K ~ Tp and T = 2.5 K 
~ 6Tp. Stated differently, we obtain an n(k) at two finite 
temperatures from the data. These n(k) are expressed 



in the form of a T = K n(k). The n(k) at the two 
temperatures are the same within experimental error. 

For a Fermi gas, the thermal broadening of n(k) is 
well known. At T ~ Tp the broadening is significant 
with n(k = 0) = 0.73, n(k F ) = 0.5, n(2k F ) = 0.23. It 
could be said that there is a rounded step of Z ~ 0.5 in 
ri(k) over a width of k ~ 2kp and n(k) develops a short 
tail. This thermal broadening is less than the resolution 
and FS broadening. The thermal broadening of n(k) for 
an interacting Fermi liquid is not known. However, in a 
strongly interacting Fermi liquid Z is already small and 
n(k) already has a long tail. Hence thermal broadening 
is expected to be relatively less important in an interact- 
ing liquid. For example, a higher temperature might be 
required before the thermal rounding (e.g. effective ther- 
mal reduction of Z) becomes significant relative to the 
'interaction' reduction. Apparently, this rounding is not 
significant or critical compared to the rounding arising 
from FS effects or the instrumental broadening within 
current error. The model parameters Z = 0.05 ±0 01 and 
77 = 0.85 ± 0.15 obtained are consistent with theory for 
T = K. Stated differently, the calculated n(k) for T = 
K reproduces our experimental data. Our observed Z 
may be somewhat smaller than the calculated value be- 
cause of resolution and thermal effects. The theory and 
experiment are therefore consistent in finding that 3 He in 
dilute 3 He- 4 He mixtures is a strongly interacting Fermi 
liquid. This strong interaction must arise at least in part 
from interaction with the 4 He component of the liquid in 
dilute mixtures. 

Mazzanti and co-workers 44 have compared their calcu- 
lated J(Q,y) for pure liquid 3 He with that observed by 
Azuah et al. at T = 1.4 K (T ~ Tp) and Q = 19.4 A" 1 
in much the same way we have done here. The input to 
their calculated J(Q,y) is n(k) for T = K calculated 
using DMC. This yields Z = 0.236 as noted above. They 
represented their calculated n(k) by an analytic function 
that had an exponential tail times algebraic terms. The 
second input is a Final State broadening function valid 
for a pure Fermi liquid. The Final State broadening func- 
tion R(Q, y) has a FWHM in y of approximately 0.7 A -1 
compared with the total FWHM of J(Q, y) of 2.0 A" 1 so 
that R{Q, y) has only a modest impact on J(Q, y). They 
find their J(Q,y) reproduces the observed J(Q,y) well 
and conclude that here is agreement between theory and 
experiment. In a similar way we have used a model n(k) 
with parameters (rather than a calculated n(k)) repre- 
sented by an analytic function with an exponential tail. 
We find a J{Q, y) (without FS effects) which reproduces 
experiment within statistical error. The model n(k) is 
consistent with n(k) calculated for mixtures^ We con- 
clude in a similar way that there is no disagreement be- 
tween theory and experiment for mixtures as well as in 
pure 3 He. 

In summary, we have measured the momentum distri- 
bution of 3 He- 4 He mixtures for 3 He concentrations x be- 
tween and 20%. We find uq increases from 7.25±0.75% 
in pure 4 He (x = 0%) to 11 ±3% for mixtures with 
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x = 15 — 20%, in agreement with theoretical calculations. 
The 4 He kinetic energy is found to decrease slightly with 
3 He concentration. The model 3 He n(k) that reproduces 
our data at x = 5% agrees with the calculated n(k) at 
x = 6.6% removing a previous apparent disagreement 
between theory and experiment based on comparing 3 He 
kinetic energies. 
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